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RESUMEN 


En este trabajo, se realiza el disefio, modelado y simulación de una bomba centrífuga, 
analizando las etapas del disefio y guiando al lector por bibliografia y metodologias sugeridas, 
con base en los parâmetros de aplicación de un proyecto real. Eltrabajo empieza con la revisión 
de literatura de disefio de bombas centrífugas y sobre antecedentes del uso de la Dinámica 
de Fluidos Computacional (CFD, por sus siglas en inglés). Posteriormente se continua con 
el disefio de los elementos que compondrán la bomba utilizando el marco teórico para 
posteriormente construir el modelo 3D en un software de Disefio Asistido por Computadora 
(CAD, por sus siglas en inglés). Al finalizar se realizan las simulaciones pertinentes dentro de 
software comercial de CFD, donde se obtendrán los resultados pare ser analizados y poder 
justificar la metodologia unificada de disefio teórico y computacional. En esta investigación se 
disefió un modelo tridimensional con base a los criterios de disefio de una bomba centrífuga 
y el uso del software Solidworks, misma que fue validadas bajo consideraciones específicas 
de las propiedades del fluido (presión y velocidad) utilizando herramientas de Mecánica de 
Fluidos Computacional con ANSYS Fluent y Solidworks FlowSimulation, obteniendo en el 
proceso una metodologia sinérgica entre la teoria con las herramientas computacionales, 
permitiendo al lector un mejor entendimiento del proceso de disefio de una bomba centrífuga. 


Resumen 


INTRODUCCIÓN 


Una gran parte de la ingeniería gira en torno a máquinas, que toman una forma 
de energia y la canalizan para ser utilizada con algún beneficio, como lo es una bomba 
centrífuga, que toma energia mecánica rotacional y la imparte a un fluido con la finalidad 
de proporcionar la energia necesaria para el flujo dentro de un sistema de bombeo, ya sea 
para distribuir un líquido, asegurar su circulación o simplemente transportarlo. 

La evolución de los sistemas de bombeo se dio debido a diferentes necesidades 
que debian ser resueltas, variantes en el tiempo al igual que las soluciones propuestas. Un 
ejemplo es el control de recursos hídrico tanto para el mantenimiento y consumo humano, 
como para el desarrollo de la agricultura y la ganadería, presentes desde las primeras 
civilizaciones. Se estima, que probablemente la primera máquina de elevación de agua 
data del 2500 a.C. construida en Egipto (Hill, 1996). El segundo sistema es el tornillo de 
Arquímedes en el siglo Illa.C., utilizado para elevar agua en pequefias alturas suministrando 
un gran caudal (Landels, 2000). Eltercer sistema utilizado, según Hill (1996), es el tympanon, 
destinado a transportar agua en alturas no muy elevadas y utilizado por los griegos mediante 
impulso por fuerza humana en el siglo Ill a.C. En cuanto a máquinas complejas para el 
transporte de agua, se tiene a la Sâgiya, construida por una cadena de cubos o jarros 
conducida por un par de ruedas dentadas para transformar el movimiento de horizontal a 
vertical (Hill, 1996), esta máquina data del siglo Ill a.C., podía funcionar con la fuerza de 
animales y era efectiva para sacar agua de los pozos a profundidades considerables. Los 
mecanismos anteriores tenían un funcionamiento meramente gravimétrico, sin embargo, 
en la Edad Antigua, también se desarrollaron bombas de desplazamiento positivo, en este 
caso, bombas de fuerza que utilizaban cilindros, pistones y válvulas, atribuídas al imperio 
romano (Landels, 2000). En cuanto al desarrollo evolutivo de los sistemas de bombeo, en la 
Edad Media no hubo un desarrollo notable en cuanto a estos sistemas (Bernal, 1979). Por 
otro lado, el desarrollo científico y tecnológico de los sistemas de bombeo siguió en oriente, 
siendo China un pionero en este campo, con la construcción de varios dispositivos entre 
los que destaca la bomba de cadena de paletas cuadradas (llamada también “columna de 
dragón” por su forma), extendiendo su uso a Corea en el siglo XV, difundiéêndose a otras 
zonas tiempo después (Needhman, 1965). 

Junto con la revolución científica en Italia en el Renaciemiento, Leonardo da 
Vinci comienza con estudio de la hidráulica de forma experimental, realizando una gran 
cantidad de artilugios, entre ellos una máquina de elevación de agua y una bomba de doble 
acción (Hill, 1996). En la primera fase del Renacimiento aparece el primer acercamiento 
a las bombas centrífugas en 1475 y se trató de una máquina para elevar barro atribuida 
a Francesco di Giorgio Martini (Moniruzzaman, 2015). Para 1588, el ingeniero italiano 
Agostino Ramelli, publica el disefo de una bomba volumétrica rotativa con mecanismo de 
paletas (Skinner, 2014 y Brashear, 2004). Los avances de estos sistemas siguen en Francia 
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y Holanda, ya que en el afio de 1600, Johannes Kepler disefió una bomba de engranes 
(Skinner, 2014). Con este invento y con las exploraciones del mundo, la tecnología empezó 
a utilizarse en diversas partes de oriente y occidente, por ejemplo, los inventos de China en 
cuanto al sistema de bombeo de agua llegaron a Holanda desde 1513 (Needham, 1965), 
donde se implementaron y mejoraron. 

En la Restauración (tercer periodo del Renacimiento), con la consolidación de las 
primeras sociedades científicas como la Royal Society de Londres y la Academia Real 
Francesa, se empezaron a dedicar a la solución de problemas técnicos de la época, 
incluidos los relacionados con la mecánica de fluidos, comenzando un significativo adelanto 
en el desarrollo de bombas, principalmente de desplazamiento positivo, como la patentada 
en 1675 por Sir Samuel Moreland, capaz de elevar grandes cantidades de agua con menos 
esfuerzo (Girdhar & Moniz, 2004). En este período, también se realizaron importantes 
avances en lo referente a bombas centrífugas, entre ellos el del físico e inventor francés 
Denis Papin, a quien se le atribuye los inicios de la investigación de bombas centrifugas, ya 
que creó un vórtice forzado en una caja circular a través de hélices (Girdhar & Moniz, 2004). 

Las bombas hidráulicas bajo la ciencia moderna, se pueden englobar en la etapa del 
capitalismo industrial, donde la ciencia y la sociedad trabajan de forma conjunta para dar 
origen al nacimiento de la civilización industrial, dando un giro importante del hombre en 
el dominio de la naturaleza, a través de la inserción de las máquinas (Bernal, 1979). Sería 
entonces hasta la Revolución Industrial en donde se consigue por primera vez introducir las 
bombas centrifugas, a partir de los hallazgos de Combs (siguiendo la teoría de Papin), en 
1838 se presenta una investigación sobre álabes curvos y el efecto curvado en los álabes, 
siendo un gran avance en el desarrollo de los rodetes para bombas (Girdhar & Moniz, 2004). 
En 1839, Andrews introdujo un disefio óptimo de la voluta uy en 1846, utilizaría un rodete 
completamente cerrado. Bajo este principio, Jhonson construyó en 1846, la primera bomba 
de tres escalamientos o de tres rodetes. Tres afios después, en 1849, James Gwynne 
construye una bomba de múltiples escalamientos, comenzando con el estudio sistémico de 
este tipo de bombas (Fraenkel, 1986). 

Si bien, el desarrollo histórico y los avances en la tecnología de bombeo han sido 
extensos y ampliamente estudiados, el disefio de estas máquinas en la actualidad es un 
proceso complejo que lleva días o incluso semanas. Este proceso se caracteriza por una 
serie de etapas: tener un problema y las condiciones a las que la máquina deberá operar, 
tomar esa información e ir construyendo un modelo que cumpla con esas características, 
el cual se someterá a pruebas, simulaciones, se realizarán las documentaciones de 
resultados, posibles fallas y mejoras, y una vez validado, se pondrá en marcha dentro el 
ambiente donde se debe desempefiar (Cola, 2017 y Covarrubias et al, 2017). 

Como parte del desarrollo y la eficiencia de tecnología de bombeo centrifugo, existen 
herramientas que hacen estos procesos más sencillos, como lo son las herramientas 
computacionales de simulación. Dichas tecnologias son esenciales en la formación del 
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ingeniero, pues ofrecen más variedad a la forma de abordar los problemas, así como 
permiten analizar y comprender conceptos que en la teoría suelen ser meramente 
abstractos, en especial si no se cuentan con modelos físicos para pruebas. 

La simulación CFD, es una herramienta para resolver numéricamente las 
ecuaciones del movimiento de los fluidos con uso de computadora (Xamán y Girón, 2015). 
Entre los trabajos más recientes de CFD aplicados al desarrollo de bombas centrífugas, 
se encuentran la optimización de métodos predictivos (Frosina et al., 2017), mejorar los 
puntos de eficiencia en el disefio (Lin et al., 2021), inestabilidad e influencias de parámetros 
geométricos (Graciano et al., 2020) y validación (Patil et al, 2020), por mencionar algunos. 

Existen actualmente una gran cantidad de softwares y herramientas para simulación 
CFD, entre ellas SolidworksO y ANSYS O. Con estas herramientas se pueden obtener datos 
cualitativos y cuantitativos importantes a la hora de validar los disefios, por lo que tener esta 
información al principio del proceso de disefio, permite tomar decisiones más integrales 
para evitar retrasos innecesarios derivados de la fabricación de múltiples prototipos físicos, 
pérdidas de presión o energia en partes importantes del modelo, y en última instancia, 
perdidas en la productividad (Solidworks Corporation, 2021). 

La metodologia empleada para el desarrollo de este trabajo, comienza explicando el 
marco teórico referente al diseio de bombas centrífugas. Posteriormente, con la información 
obtenida en el marco teórico, se procederá a realizar los cálculos de los elementos que 
integran la bomba, tomando como base los parâmetros y variables de entrada y salida con 
referencia a un problema de aplicación. Mediante el uso de herramientas de modelado CAD 
como es SolidworksQ, se desarrollarán los modelos 3D de los diferentes componentes 
que conforman la bomba centrífuga. Una vez desarrollados todos los componentes, se 
realizará el ensamblaje completo de la bomba centrífuga. Para finalizar se usará ANSYSO, 
un software de ingeniería avanzada para la simulación fluidodinámica del ensamble de la 
bomba centrifuga. 

El objetivo de este trabajo, es proporcionar una metodologia para el disefio de una 
bomba. Se analizará la bibliografia recomendada y mediante un caso de estudio se aplicará 
lo investigado y las tecnologias computacionales propuestas. Finalmente se contrastarán 
los resultados para determinar la viabilidad de utilizar métodos combinados de la teoría y 
cálculos simples con los softwares de simulación comerciales. 
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ANTECEDENTES 


Para satisfacer las necesidades de la industria, los servicios y muchas actividades 
humanas, se necesitan máquinas que realicen distintas operaciones, una de ellas es la 
bomba centrífuga, que entra en la clasificación de máquinas de fluidos y posteriormente 
en las hidráulicas. Considerando a Garcia (2020), cada máquina es un sistema de 
componentes mecánicos que intercambia energia, pasándola de una forma a otra. Las 
máquinas de fluidos realizan esta operación mediante líquidos y gases. Por su parte las 
máquinas hidráulicas son una variedad de máquina de fluidos en las que los efectos de 
compresibilidad del fluido son despreciables, es decir trabaja con líquidos de volumen 
constante. 

Asu vez y concordando con Fontalvo Conrado et al. (2020), las máquinas hidráulicas 
pueden ser divididas en dos tipos esenciales, ya sea las turbinas que obtienen o extraen 
energia del fluido en cuestión con la finalidad de convertirla en energía útil, como puede 
ser los generadores eléctricos; o las bombas que hacen exactamente lo opuesto. Toman 
energia de una fuente y la convierten para impartírsela al líquido de trabajo y producir un 
flujo másico aumentando la velocidad de las partículas del fluido. La bomba es la parte 
fundamental de un sistema de bombeo para distribuir la sustancia. 

La bomba centrífuga es la más común en todas las aplicaciones La constante 
necesidad de sistemas de bombeo y las múltiples aplicaciones que tienen en la industria 
y sectores residenciales. Por tanto, es importante conocer y entender su funcionamiento 
y los principales factores que intervienen en su desempefo, con el fin de obtener disefios 
óptimos desde un punto de vista energético. 

Para disefiar una bomba, se requiere amplio conocimiento del flujo interno, así como 
la forma en que influyen las condiciones de operación. Con suficiente conocimiento se 
puede realizar un disefio basândose en la teoría y cálculos metódicos que hacen el proceso 
rápido, pero algo alejado de la realidad. En contraste se utilizan los modelos y las pruebas 
experimentales, no obstante, el proceso es tedioso y toma grandes porciones de tiempo y 
recursos al tener que fabricar varias piezas y realizar muchas pruebas. 

Según lo dicho por Shah et al. (2013) en los afios recientes, la Dinámica de Fluidos 
Computacional (CFD) ha dominado ambas partes, pues ofrece resultados muy cercanos a 
la realidad en intervalos de tiempo considerablemente rápidos. Hace casi 40 afios comenzó 
a surgir esta rama de dinâmica de fluidos. Es de esperarse que su desarrollo fue iniciado 
por la introducción de teoría del elemento finito en un entorno computacional. 

García (2020) menciona que en general el análisis por elemento finito es dividir el 
sistema en elementos más pequefios interconectados por nodos a lo largo de una malla. 
Naturalmente entre menor sea el tamafio de los elementos, más real será el cálculo y, 
por lo tanto, más preciso; sin olvidar que también habrá un mayor número de elementos, 


haciendo el proceso más largo y complejo. De esta forma se debe buscar un balance entre 
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exactitud y complejidad. 

El análisis computacional de las bombas centrífugas se hace generalmente mediante 
software comercial como lo es Fluenty ANSYS-CFX. Pinto et al. (2017), establece que se 
deben considerar una gran cantidad de parâmetros y como se relacionan entre sí. 

De acuerdo a lo comentado por Spence & Amaral-Teixeira (2009) y Medvitz et al. 
(2002) el desempefio de una bomba está fuertemente ligado a su geometria y un sinfín de 
condiciones, como lo son la distribución de los álabes, el espacio donde comienza la voluta 
y donde termina, las paredes internas de la bomba y sus obstrucciones, y el espesor de los 
álabes. Otro componente importante es el de la cavitación, que ha sido investigado, donde 
se pueden analizar los efectos que este fenómeno conlleva en una máquina de bombeo. 

Mediante el uso de CFD podemos obtener un método de tres etapas para analizar 
el desempefio de una bomba centrífuga dentro de las condiciones deseadas. La primera 
etapa es generar una malla adecuada para el modelo, que incluye los componentes sólidos 
y líquidos. Posteriormente en la segunda fase, se deben descubrir las condiciones de 
frontera e iniciales con las que se realizará el estudio, esto dictará el comportamiento del 
fluido. Y la etapa final es realizar distintas simulaciones para obtener la mayor información 
posible y poder destacar los puntos clave de disefio. 

Según Yang et al. (2011), utilizando modelos matemáticos de volumen finito se 
puede conjugar la ecuación de transporte para disefiar una voluta con la máxima eficiencia 
esperada. No obstante, este método demostró que eficiencias de mayor valor no pueden 
ser simplemente encontradas con la optimización de los cálculos de disefo. 

Pinto et al. (2017) menciona que si bien la dinámica de fluidos computacional promete 
ser una herramienta compleja pero muy útil, conlleva sus desventajas, la principal siendo 
su alta demanda de tiempos y recursos de cómputo. Obtener simulaciones avanzadas 
requieren equipos de cómputo más robustos; a su vez los tiempos de cálculo de algunos 
modelos pueden tomar días y esto gradualmente impacta en costos. 


11 JUSTIFICACIÓN 


Los programas de Disefio asistido por computadora (CAD, por sus siglas en inglés) 
son herramientas importantes que los estudiantes deben dominar para aumentar el 
potencial de disefio. Al implementar una investigación sobre el contraste y las similitudes 
entre la parte teórica y los simuladores de flujo, se espera lograr una metodologia confiable 
de validación para ser implementada en la parte académica del Instituto Tecnológico de 
Hermosillo, especificamente en las asignaturas referentes a la mecánica de fluidos. 
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21 OBJETIVOS 


2.1 Objetivo general 


Disefiar un modelo tridimensional con base a los criterios de diseho de una bomba 
centrífuga y el uso del software Solidworks, misma que será validadas bajo consideraciones 
específicas de las propiedades del fluido (presión y velocidad) utilizando herramientas de 
Mecánica de Fluidos Computacional con ANSYS Fluent y Solidworks FlowSimulation, 
obteniendo en el proceso una metodologia que concilie la teoria con las herramientas 
computacionales y permita al lector un mejor entendimiento del proceso de disefio de una 
bomba centrífuga. 


2.2 Objetivos específicos 


*— Investigar el marco teórico referente a los criterios de cálculos para bombas 
centrifugas. 


. Definir los criterios y características de operación de la bomba. 


* Realizar los cálculos de disefo para los elementos esenciales que conforman 
la bomba. 


* | Modelar los elementos de la bomba en software de disefio. 
* | Ensamblar los elementos para obtener la bomba en un disefio tridimensional. 
* Realizar los estudios de CFD para la bomba modelada en 3D. 


*— Analizar los estudios de CFD para posibles cambios de optimización en el di- 
serio. 


* Validación del diseio a partir de resultados teóricos y de CFD. 


2.3 Alcances y limitaciones 


La investigación realizada abarcará el uso de dos âmbitos principales, el primero 
son los cálculos de disefio teóricos obtenidos de la bibliografia recomendada de bombas 
centrífugas y segundo el uso de software de disefio asistido por computadora, así como 
dinámica de fluidos computacionales. No se realizarán modelos físicos ni se harán pruebas 
reales, la experimentación será meramente computacional. 

Se propone una situación teórica de ejemplo sencilla donde se recoge la información 
principal para el disefio de la bomba y se guiará por la metodologia de disefio basado en la 
literatura recabada en el marco teórico. Los cálculos realizados serán a nivel universitario 
con ecuaciones que se ven generalmente en cursos de dinámica de fluidos, sistemas 
hidráulicos y máquinas de fluidos incompresibles, esto con la finalidad de servir como guía 
accesible a estudiantes de nivel licenciatura. 

El disefio será realizado mediante la plataforma de Solidworks pues ofrece una 
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gran variedad de herramientas útiles para realizar los modelos virtuales manteniendo 
una interfaz amigable para el usuario. Las simulaciones se calcularán dentro de dos 
softwares de dinâmica de fluidos, SolidworksO Fluid Simultaiony ANSY SO Fluent usando 
solamente las herramientas básicas que ofrecen ambas plataformas y con configuraciones 
recomendadas por las mismas. 


2.4 Hipótesis 


Los cálculos teóricos permiten una aproximación sencilla pero no tan fiel de la 
realidad, pero al ser utilizados como base de una simulación se obtendrán resultados más 
reales y convincentes para la validación de una bomba centrífuga. 


2.5 Planteamiento del problema 


El disefiar una bomba centrífuga puede ser algo complejo y tedioso si no se cuenta 
con una metodologia y herramientas adecuadas. El reto en cuestión es en consultar 
bibliografia reciente y confiable, así como determinar cuál es el mejor camino al momento 
de realizar la comprobación del disefio. También resulta complejo el adentrarse en la 
complejidad de los sistemas de dinámica de fluidos computacional sin haber tenido un 
acercamiento previo y sin conocer los conceptos principales. 

La problemática principal es el poco acceso de un caso práctico que guíe al lector 
por el proceso de disefio y sus etapas, que también se demuestre que los cálculos teóricos 
y la simulación no son dos opciones diferentes por las cuales uno puede optar al momento 
de resolver un problema, sino herramientas que en conjunto permiten la solución de los 
problemas de una forma más eficaz y con mejores resultados. 

Es de esto que surgen las dos principales preguntas de investigación para el 
proyecto: 

* | é Son coincidentes los resultados que se obtienen de forma teórica con respecto 


a los que se obtienen mediante simulaciones de Dinámica de Fluidos Compu- 
tacional? 


* é Cómo es una metodologia de disefio que utilice la facilidad de la teoría y tenga 
los beneficios de las herramientas computacionales? 
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MARCO TEÓRICO 


El principio de funcionamiento de una bomba centrífuga, como su nombre lo 
indica, se basa en transmitir energia desde una cámara con un componente rotacional, 
que transmite la energia hacia el fluido, provocando su rotación y expulsándolo mediante 
la fuerza centrífuga resultante, posteriormente el fluido se recoge por una espiral y es 


canalizado por una câmara convirtiendo su energia cinética en energia de presión. 


11 ELEMENTOS CONSTITUTIVOS DE LA BOMBA CENTRÍFUGA 


La bomba centrífuga tiene dos componentes principales que realizan las dos 
funciones principales que son impartir energía al fluido y convertir la presión generada a 
velocidad, estas partes son el impulsor y la voluta respectivamente. 


1.1 Impulsor 


La parte central de la bomba centrífuga, un elemento rotativo que succiona el líquido 
por su centro y le imparte velocidad al fluido al guiarlo por un conjunto de álabes. Este 
componente tiene muchas variaciones y es el punto central del disefio de una bomba. 

Se pueden distinguir en distintos dipos de impulsores: 

* En un impulsor abierto, los álabes desnudos van unidos únicamente al eje de 


giro y se mueven entre dos paredes laterales fijas pertenecientes a la carcasa 
de la bomba. 


* Los impulsores con una sola pared lateral, que siempre es la posterior, se em- 
plean con cierta frecuencia, destacando las bombas de flujo mixto y todas las 
axiales. 


* | Los impulsores cerrados tienen los álabes colocados entre dos paredes late- 
rales, anterior o de aspiración y posterior 


* | Los impulsores de doble aspiración llevan aros de cierre en los dos oídos. Se 
pueden considerar como dos impulsores de aspiración simple, opuestos y en 
paralelo. 


* Los impulsores de aspiración simple cuando están provistos en la parte poste- 
rior de câmara de equilibrado del empuje hidráulico axial en comunicación con 
la aspiración a través de los agujeros de equilibrio, sólo tienen aros a ambos 
lados 

De acuerdo a Fernández Díez (2012), los álabes del rodete someten a las partículas 
de líquido a un movimiento de rotación muy rápido, siendo proyectadas hacia el exterior por 
la fuerza centrífuga, de forma que abandonan el rodete hacia la voluta a gran velocidad, 
aumentando su presión en el impulsor según la distancia al eje. La elevación del líquido 
se produce por la reacción entre éste y el rodete sometido al movimiento de rotación; en 
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la voluta se transforma parte de la energia dinámica adquirida en el rodete, en energia 
de presión, siendo lanzados los filetes líquidos contra las paredes del cuerpo de bomba y 
evacuados por la tubería de impulsión 


Cerrado De doble aspiración Semiabierto Abierto 


Ilustración 1. Distintos tipos de Impulsores 


1.2 Carcasa 


Es el componente que permite convertir la energia cinética del fluido proporcionada 
por el impulsor, en energia de presión. Esto se logra con una reducción gradual de la 
velocidad del fluido aumentando el área de flujo. La carcasa o voluta está dispuesta 
en forma de caracol, de tal forma que la separación es mínima en la parte superior; la 
separación aumenta hasta que las partículas del fluido se encuentran frente a la salida 
donde está la descarga. 

De acuerdo a Parra & Robles (2016) el sistema difusor de una bomba consta de tres 
elementos: 

* | Corona directriz. Tiene dos finalidades, sirve para guiar al agua hacia la caja 
espiral o voluta de la forma hidráulicamente más eficiente, también para dismi- 
nuir la energia cinética que tiene el fluido en la salida del rodete y convertirla 
en presión 


* | Cajaespiral o voluta. Tiene una doble función; por un lado, recoge el fluido de la 
periferia del rotor o del difusor y lo encauza hacia la conducción de salida, y por 
otro, actúa también como un difusor, puesto que tanto su radio como su sección 
creciente disminuyen la velocidad y aumentan la presión. 


* | Cono difusor. Es la parte final, una terminación cónica que se conecta mediante 
una brida a la tubería de descarga. 


No siempre se encuentran los tres elementos, como en la figura podemos ver 


distintos difusores con diferentes características. 
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llustración 2. Difusores. A) Con caja espiral. B) Con caja espiral y cono difusor. C) Con caja espiral, 
cono difusor y corona directriz. 


1.3 Otros elementos 


Existen más elementos que constituyen una bomba centrífuga, entre ellos podemos 
mencionar las bridas de succión y descarga que se encuentran a la entrada y la salida de 
la bomba. La brida de succión conecta el interior de la bomba con la tubería por donde el 
fluido es succionado hacia el impulsor. En contraste la brida de descarga es la terminación 
de la voluta y conecta está a la tubería de descarga, donde el fluido emerge de la bomba 
hacia el sistema de tuberías. 

Gómez Osorio & Zarate Pernett (2004) comentan sobre la flecha que es otro elemento 
de la bomba que tiene como función el transmitir los torques que se presentar al arrancar 
y durante la operación, a su vez soportando el impulsor y el resto de los componentes 
giratorios. Las flechas deben soportar el esfuerzo que se aplica al arrancar súbitamente 
una bomba. 


21 FLUJO DENTRO DE LA BOMBA CENTRÍFUGA 


Dentro de una bomba centrífuga la energia cinética del rodete se transmite al fluido, 
generando un cambio en la velocidad de este. Para realizar cálculos del fenómeno, usamos 
la ecuación de flujo de Euler, utilizando diagramas vectoriales con el fin de visualizar el 
comportamiento. 


2.1 Triángulos de velocidades 


Las velocidades dentro del impulsor se pueden demostrar mediante triângulos de 
velocidades, que relacionan tres velocidades, la velocidad absoluta del fluido (c, a entrada 
y c, a salida) y sus dos componentes, una velocidad tangencial del impulsor (u, y u,) y una 
velocidad relativa al movimiento del impulsor (w,y w,). También se incluyen los ângulos a 
la entrada y salida de los álabes (B, y B,) que son los ángulos entre la velocidad absoluta y 
la tangencial; así como los ángulos entre la velocidad relativa y la tangencial (a, y a,) y los 
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ângulos entre la velocidad absoluta y la relativa (y,y y,). 

Es importante mencionar que se consideran las proyecciones de las velocidades 
respecto a la tangencial, es decir la velocidad absoluta y la velocidad relativa proyectadas 
a la tangencial (c,y w,). Finalmente, la velocidad meridional, que es una proyección de 
la velocidad absoluta perpendicularmente a la tangencial (c ), esta se le conoce como la 
velocidad meridional, es perpendicular a la circunferencia del impulsor. 

Si bien estas velocidades no se conocen exactamente dentro del impulsor, se 


pueden hacer aproximaciones utilizando las siguientes expresiones: 


2nnr 
u= 1 
60 
Qa 
= 2 
m  2mrb 


Utilizando: 

n = la velocidad de giro del impulsor en rpm, 
r = radio respecto al centro del impulsor, 

b = el ancho por donde pasa el fluido, 

Q,=y elflujo de trabajo del impulsor. 


Ilustración 3 Diagrama vectorial de velocidades a la entrada y salida del impulsor. 


2.2 La ecuación de Euler 


La ecuación de Euler expresa que la potencia transmitida al fluido por unidad masa 
es: 
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Donde: 

Y, es la potencia trasmitida por unidad de masa; 

P.. potencia teórica que transmite el impulsor al fluido; 
Q.. el flujo que pasa por el impulsor; 

p, densidad del fluido; 

c 


C., Son las proyecciones de la velocidad del fluido respecto a la velocidad 


2u? “tu? 


tangencial a la salida y entrada, respectivamente. 

No obstante, es más común referirse a la altura que proporciona la bomba, es decir, 
energia por unidad de peso de fluido (H), que se expresa de la forma: 
LO — Coulz — Gu 
9 9 


Si analizamos esta expresión, se puede notar que se puede maximizar la altura 


Ha = 4 


cuando se obtiene c, =0, es decir cuando el ángulo a=90º, es decir la entrada del fluido es 
perpendicular. Para esto se implementan álabes directrices en la tubería de aspiración y así 
evitar que el fluido Ilegue con pre-rotacion. 

Es importante mencionar que en la práctica no se conocen las verdaderas 
velocidades ni sus direcciones. Los triângulos de velocidad se dibujan sobre los ángulos 
de los álabes y sobre eso y demás suposiciones se va desarrollando hasta completar el 
triângulo y obtener un cálculo de la altura de Euler. 


2.3 Eficiencia de la bomba 


Es importante mencionar que no toda la potencia que proporciona el impelente se 
convierte en potencia útil al fluido. Todas las pérdidas que se originan entre los puntos 
donde se mide la presión de succión y descarga, constituyen las pérdidas hidráulicas. 

Según Marchegiani (2004) estas incluyen pérdidas por fricción a lo largo de la 
trayectoria del líquido desde la brida de succión hasta la de descarga; pérdidas debidas a 
cambio brusco, tanto en área como en dirección de flujo; y todas las pérdidas debidas a 
remolinos, cualquiera que sea su causa. Por lo tanto, se utiliza una eficiencia para ajustar 
a estas pérdidas, la eficiencia hidráulica n, 


[o 
pgla 


Además de las pérdidas hidráulicas, Marchegiani (2004) dice que existen pérdidas 


H=mHa=m = Na(Cullz E Cruls) 5 


volumétricas, debido a las fugas que existen en los espacios entre partes rotatorias y 
estacionarias de las bombas. El caudal en la descarga de la bomba es menor que en la 
succión y también, es menor que el caudal que pasa por el impulsor. El cociente de los dos 
caudales de descarga y el que pasa por el impulsor se Ilama rendimiento volumétrico. 
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Da 0a 
Q Qa + Qi 


Con n, representando la eficiencia volumétrica y Q, el volumen de las fugas. 


No 


Finalmente se tiene la eficiencia mecánica que es la relación que existe entre la 
potencia que realiza la bomba al suministrado por el eje. Este se debe a componentes 


como cojinetes, sellos y la misma fricción de los componentes. 


De esta forma, utilizando todos los coeficientes de las eficiencias obtenemos la 
potencia final también Ilamada potencia al freno y lo relacionamos con la teoría de Euler. 
— P9QHa — pgQH 

MaNvAm Ne 


P 8 


Siendo n=n,n,n,, es decir, la eficiencia total de la bomba. 
Si bien las eficiencias sirven meramente para acercar los cálculos más a lo real, no 
son la mejor opción para calcular las pérdidas reales de una bomba, por lo tanto, no son 
tan viables para el disefio pues se basan en suposiciones. No obstante, existe bibliografia 
aceptada respecto a los valores comunes que estas pueden presentar. 
Existen diferentes motivos por los cuales no se presentan la altura de Euler, entre 
las que se pueden resumir: 
1. La velocidad relativa del fluido en la parte posterior es mayor que en la parte 
frontal del impulsor. Esto se debe a las distribuciones de presión dentro del mismo. 


2. Existe un fenómeno conocido como circulación relativa, que se ocasiona por 
efectos de inercia en el agua por el giro de las regiones cercanas, resultando en un 
componente opuesto a la velocidad tangencial y por lo tanto reduciendo el valor de 
la altura útil. 


2.4 Corrección por el número de álabes 


Las ecuaciones antes mencionadas funcionan dentro de rodetes teóricos con un 
número infinito de álabes, para adecuar esta teoría a la realidad, se utilizan correcciones 
que ajustan los datos a la desviación por los mismos álabes. Se utilizará la corrección de 
Pfleiderer, que estipula lo que la altura con álabes infinitos (H.) y la altura H, con Z álabes 
se relacionan de la forma: 

Ho 


ci NI — (11/15)*] 


Siendo el coeficiente de Píleiderer y=(0.55/0.65)+0.6senB, 
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2.5 Coeficiente de obstrucción 


Otra corrección importante a la teoría ideal de las bombas centrífugas es el uso de 
un coeficiente 1 que resulta del espesor de los álabes. Esto se debe a que el espesor de 
los álabes produce una reducción en el área a la entrada y la salida del impulsor. Dicho 
coeficiente se puede obtener con la ecuación: 


E TD — Z(e/senB) 


10 
TD 


Que siempre resultará ser menor a 1, de esta forma podemos corregir el flujo que 
pasa por el impulsor, ya sea en la entrada o salida, solo multiplicando el área por t. Cabe 
destacar que en algunos casos los álabes son afilados, es decir tienen un espesor mínimo, 
en estos casos el coeficiente toma el valor de 1. 


31 CLASIFICACIÓN DE LAS BOMBAS CENTRÍFUGAS 


En orden de obtener una clasificación es necesario acudir a un sistema que 
permita agrupar a bombas de la misma naturaleza, de ahí la introducción de un coeficiente 
denominado velocidad específica o, en bibliografias más antiguas, número específico 
de revoluciones. Para dos bombas con una velocidad específica similar, se dice que son 
semejantes, puesto que se espera un comportamiento igual en respuesta a las condiciones 
de operación. La velocidad específica O se define como: 

nQ!2 
És = cgHyIA di 

Es decir, es una función de las revoluciones (O) la altura útil y el caudal en el punto 
de máxima eficiencia de la bomba. Tiene valores adimensionales y ronda entre 0.2 y 2.0 
para las bombas centrífugas. 

Barreto (2004b). comenta en la velocidad específica que es de amplio uso en el disefio 
de las turbomáquinas. Permite establecer un prototipo a partir de un modelo con formas, 
proporciones geométricas y velocidades conocidas. Originalmente fue implementado con 
éxito en el disefio y construcción de turbinas y su uso se extendió posteriormente a bombas, 
ventiladores y compresores. 

Otro número adimensional común para el estudio y disefio de bombas centrífugas 
es el diámetro específico A, que se puede calcular de la forma: 


— D(gh)4* 


Este par de números adimensionales funcionan muy bien para caracterizar bombas 
centrífugas, no obstante, resultan útiles para disefiar bombas de alta eficiencia, utilizado 


Marco Teórico 


15 


en distintas bibliografías. El método más común para el disefio es el diagrama de Cordier. 


50 
40 
30 


Ilustración 4. Diagrama de Cordier. 


3.1 Tipos de bombas centrífugas 
Las bombas centrífugas pueden ser clasificadas en: 
* Centrífuga o de flujo axial 
*— Flujo mixto 
* Flujo radial 


En las bombas centrífugas radiales la corriente líquida se verifica en planos radiales, 
en las axiales en superficies cilíndricas alrededor del eje de rotación y en las diagonales se 
verifica radial y axialmente, denominâándose también de flujo mixto. 

La forma de los álabes en los impulsores de flujo radial es, en general, curvada hacia 
atrás con respecto al sentido de giro, B2 < 90º, y con superficies de simple curvatura, siendo 
la generatriz paralela al eje de rotación; en los impulsores helicoidales, los álabes son de 
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doble curvatura y en los axiales tienen, además, un determinado perfil aerodinámico. 


41 CAVITACIÓN 


Marchegiani (2004) comenta que el fenómeno de la cavitación es muy importante en 
las bombas ya que será el factor determinante del funcionamiento del sistema de bombeo. 
Es sabido que cuanto mayor sea la velocidad de rotación, mayores han de ser las pérdidas 
de carga que se originen a la entrada del impulsor. Estas pérdidas a la entrada de la bomba 
pueden resultar peligrosas. La presión en cualquier punto de la línea de succión dentro de 
un sistema de bombeo nunca debe reducirse a la presión de vapor del líquido. 

Cuando la presión del fluido se reduce a su presión de vapor, este se vaporiza en 
forma de burbujas, las cuales son arrastradas junto con el líquido hasta una región donde 
se alcanza una presión más elevada y allí desaparecen; a este fenómeno se le conoce 
como cavitación y tiene varias consecuencias detallas a continuación. 

Las burbujas de vapor que disminuyen el espacio utilizable para el paso del líquido, 
se perturba la continuidad del flujo debido al desprendimiento de gases y vapores disueltos, 
impactando en el rendimiento de la bomba; en su recorrido darfian los con- ductos de paso 
del líquido en el tubo de aspiración y Illegan a una zona en el rodete, de presión superior a 
la presión de vapor, en la que, instantâneamente, toda la fase de vapor pasa a líquido, de 
forma que el volumen de las burbujas pasa a ser ocupado por el líquido, en forma violenta, 
que se acompafia de ruidos y vibraciones, lo cual se traduce en un golpeteo sobre los 
álabes, que se transmite al eje, cojinetes, cierres mecânicos, etc. 

Fernández Díez (2012) menciona que si la bomba funciona en estas condiciones 
durante cierto tiempo se puede dafiar; la intensidad del golpeteo a medida que disminuye 
la presión absoluta a la entrada del rodete. 


51 NPSH 


Mott Robert L. (2006) menciona en su libro de dinámica de fluidos que los fabricantes 
y disefiadores de bombas prueban las bombas para determinar el nivel de succión 
requerido, es decir la presión necesaria para evitar la cavitación. Esto se reporta como dato 
de la bomba bajo el nombre de Altura Neta de Succión Positiva Requerida (Required Net 
Positive Suction Head, NSPH,) Básicamente este valor engloba la presión de vaporización 
del fluido y las pérdidas que el equipo tiene bajo cada punto de operación. Al momento 
de seleccionar una bomba resulta esencial este valor para no sufrir de los impactos de la 
cavitación al no permitir que la presión de succión se reduzca a dicho valor. 

En contraste, para la aplicación de una bomba es el valor de la Altura Neta Positiva 
de Aspiración (Available Net Positive Suction Head, NSPH,) el cual se determina como la 
energia disponible que existe desde el pozo de succión hasta la entrada de la bomba. Para 
definir esto es necesario determinar la energia bruta disponible a la entrada, utilizando la 
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ecuación de Bernoulli se puede calcular. Es responsabilidad del disefiador del sistema de 
bombeo que la altura neta disponible sea bastante mayor que la requerida. El Instituto 
Nacional Americano de Estándares (ANSI, por sus siglas en inglés) recomienda como 
margen mínimo 10% del valor del NSPH, sobre el NSPH,. En ciertas aplicaciones se 
esperan márgenes más elevados, como del 100%. Es claro que entre mayor sea el exceso, 
menor posibilidad existe de enfrentar los problemas ocasionados con la cavitación, a lo cual 
también contribuye altamente la documentación y el criterio de selección del dispositivo de 
bombeo. 


61 MATERIALES PARA BOMBAS CENTRÍFUGAS 


Existen varias consideraciones para realizarse al seleccionar un material para 
bombas centrífugas. Podemos mencionar la cavitación, que es un fenómeno donde las 
pérdidas de presión ocasionan que el fluido pase a un estado líquido y genere burbujas 
explosivas que desgastan el material. También, debemos hablar sobre estabilidad química 
y corrosión, además de considerar impurezas en el fluido que puedan ocasionar abrasión 
dentro de la bomba. 

Según Barreto (2004a) la cavitación puede evitarse asegurândose que exista 
suficiente presión a la succión de la bomba, no obstante, la corrosión y abrasión son factores 
muy importantes para considerar. Para los rodetes sometidos a las condiciones agresivas 
del líquido tiene importancia la resistencia a la corrosión, lo que, en última instancia, define 
la vida útil de la máquina en ausencia de cavitación 

Los aceros al carbono son los más utilizados para fabricación de impulsores, debido 
a sus propiedades mecánicas y buena resistencia a la corrosión. También los aceros 
inoxidables son apropiados para aplicaciones en medios ácidos y corrosivos. Otro material 
común son las aleaciones de bronce, que tienen mayor estabilidad química, permiten 
mayor durabilidad. 

Desde los afios cincuenta transcurre la introducción de los plásticos en la fabricación 
de diferentes elementos de las bombas centrífugas. Este proceso ha avanzado de forma 
relativamente lenta, debido a que las propiedades mecânicas y físicas de estos materiales 
no siempre responden a las exigencias de distintos elementos de máquinas y, en gran 
medida, dependen del tiempo de acción de las cargas. Por otro lado, estos materiales están 
limitados a temperaturas relativamente bajas y, además, la tecnología de fabricación tiene 
sus propias exigencias que, en ocasiones, obliga a cambiar las formas más adecua- das. 
En (Barreto, 2004a) se establece que las máquinas con velocidades de rotación más bajas 
y, en consecuencia, menores velocidades específicas, se ven compensadas en parte, por 
mejores condiciones de aspiración. 
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6.2 Simulación de bombas centrífugas 


La dinámica computacional de fluidos es unarama de la ciencia e ingeniería que estudia 
el comportamiento de los fluidos haciendo uso de herramientas computacionales Dicha 
rama comprende el estudio del comportamiento de fluidos, el modelamiento matemático y 
discretización de estos modelos y la implementación en lenguaje computacional para el uso 
de dichas herramientas. El CFD permite realizar modelamiento y optimización de distintos 
procesos en los cuales el comportamiento del fluido es crítico. 

Pineda Pérez (2014) afirma que las ecuaciones que gobiernan CFD son la ecuación 
de continuidad y las ecuaciones de conservación de momento y energia. Estas ecuaciones 
se derivan de adoptar las leyes físicas de conservación de masa, segunda ley de Newton 
y primera ley de la termodinámica. Otras ecuaciones adicionales utilizadas en CFD son las 
16 ecuaciones para modelar el flujo turbulento, el cual se encuentra presente en gran parte 
de los problemas de ingeniería. 

La gran mayoría de los programas comerciales utilizan el método de volúmenes 
finitos para resolver las ecuaciones mediante la Mecánica de Fluidos, en general se plantea 
lo siguiente: 


* El dominio se divide en un número muy grande pero finito de volúmenes de 
control, que no son más que celdas con volumen determinado. 


«| En cada celda se resuelven las ecuaciones generales de conservación de 
masa, cantidad de movimiento, energia y otras. 
De esta forma, Ballesteros (2003) comenta que las plataformas comerciales permiten 
el acceso sencillo a estos algoritmos mediante interfaces amigables para el modelado de 
geometrias y la introducción de parámetros de solución 


6.1 Preparando un problema en dinámica de fluidos computacional 


Se puede resumir en 5 etapas el proceso de simulación mediante CFD, como 
aparece en el libro Computacional Fluid Dynamics: A Practical Approach. De Tu, Jiyuan 
(2018) 

El primer paso en cualquiera análisis de CFD es generar la geometria de la región 
de flujo, a esto se le conoce como dominio computacional. Dentro de este dominio se 
realizarán todos los cálculos respecto al fluido, sus propiedades, su dinâmica y el medio 
que lo rodea. Los dominios pueden ser bidimensionales o tridimensionales, dependiendo 
del tipo de análisis que se desee realizar. 

La siguiente etapa del proceso de análisis es el generar un mallado adecuado para 
el dominio, podría considerarse la parte más importante de la configuración previa a la 
simulación. La premisa es que el dominio previamente establecido debe ser subdividido 
en un número entero de geometrias más pequefias, no encimadas y relacionadas entre 


sí, las cuales serán resueltas y conllevarán a una solución final del dominio completo. 
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Este proceso da como resultado una malla, que es un conjunto de celdas cuyas aristas se 
encuentran en nodos. 

La precisión de la solución está fuertemente relacionada con la cantidad de celdas en 
las que se divide el dominio. Aunque esto sea relativamente cierto, hay otras propiedades 
que influyen en la precisión, como lo son el tipo de mallado, la exactitud de los modelos 
matemáticos y el método de solución seleccionado. 

De acuerdo con Tu (20018) la tercera parte consiste en determinar los procesos 
físicos relacionados con el análisis. Primero se debe hablar si el sistema está en equilibrio 
o si entramos en un estado transitorio donde este evoluciona respecto al tiempo. A partir de 
esto debemos considerar las propiedades del fluido, si es viscoso o no; también que tipo de 
flujo es, ya sea laminar o turbulento; si hablamos de un material compresible o incompresible. 
Además, se pueden agregar otros fenómenos como elementos electromagnéticos o de 
transferencia de calor. 


6.2 Condiciones de frontera 


El cuarto paso hace referencia a las condiciones que se imponen al sistema para 
el problema que se debe resolver. Se trata de determinar las condiciones que dentro de la 
simulación son permisibles y están disponibles como suposiciones para el dominio de flujo. 
Pueden ser tan simples como entradas y salidas de flujo másico o volumétrico; así como 
puede haber entradas de presión, presiones de ambiente; en contraste tenemos paredes ya 
sean en movimiento o estáticas que pueden ser obstáculos o impartir movimiento al fluido. 
También existen fronteras abiertas donde las propiedades del fluido se liberan al entorno. 

Es de gran importancia determinar las condiciones de frontera adecuadas para 
el modelo, puesto que más precisas sean, el cálculo convergerá de forma más rápida 
ofreciendo mejores resultados en tiempos menores. Esto se puede comprender gracias a la 
idea de la información de entrada y de salida; entre más información se tenga del fenómeno 
estudiado, más rápido se pueden encontrar las soluciones de este. 


6.3 Solución numérica 


Para utilizar adecuadamente un software de CFD es importante conocer el código 
numérico que este utiliza para realizar los cálculos y obtener los resultados. Esto se refiere 
al solver que es un algoritmo iterativo que provee soluciones, las revisa y corrige hasta 
converger a un estado que se considera la solución correcta. 

Se puede analizar un solver como tres partes, inicialización, control de solución, 
monitoreo de solución, cálculos CFD y verificación de convergencia. En este orden lo 
primero es proponer una solución inicial y determinar los parâmetros y variables para los 
cuales se resuelve, así como los criterios para una solución aceptable, se proceden a 
realizar los cálculos y si estos no son satisfactorios respecto al a convergencia, se repite el 
proceso con nuevos datos a partir de los anteriores. 
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Início 


Inicialización 


Control de 
solución 


Monitoreo de 
solución 


Cálculo de CFD 


Modificar la 
solución 


Verificar 
convergencia 


Resultados de 
simulación 


Final 


Ilustración 5. Diagrama de flujo del algoritmo del Solver 


Si las condiciones iniciales son cercanas a la respuesta, la solución será rápida; 
caso contrario si las variables iniciales son distantes de lo real, el solvertomará más tiempo. 
Los métodos matemáticos para solucionar los modelos son complejos y consideran todas 
las variables del fluido, pero se enfocan en ciertos tipos de cálculo para especializarse en 
sistemas diferentes. Los programas de CFD tienden a contar con algoritmos optimizados y 
especializados para distintas clases de simulaciones. 

El quinto paso se presenta al realizar los procesos de solución. Cuando el solver 
calcula una iteración tiene como objetivo encontrar el valor de una variable que será 
monitoreada además de resolver las variables de operación de cada celda. A su vez calcula 
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errores que determinan la cercanía hacia la solución real. Cuando después de varias 
iteraciones, los valores tienden a converger a un valor en concreto y es entonces cuando 
se puede decir que se encontró una solución final, siempre y cuando se encuentre dentro 
de los rangos de error aceptados. 


6.4 Resultados y post-proceso 


Es común asociar los programas de CFD con la producción de imágenes vívidas de 
resultados de las variables dentro del sistema simulado, que si bien son muy promocionales, 
resultan altamente útiles al analizar el modelo resuelto. Ballesteros (2003) explica que 
por esto los modelos comerciales contienen métodos gráficos avanzados y altamente 
concentrado en los resultados de cada sección mediante recursos visuales Ilamativos pero 
fáciles de interpretar. 

También se pueden reportar resultados de otras formas, como lo son tablas 
bidimensionales de interrelación entre dos variables, gráficos de vectores dentro de un 
plano, gráficas de gradientes, líneas de flujo. También, dado que son entornos multimedia, 
también pueden introducirse animaciones, así como se pueden manipular las vistas del 


modelo. 


Configuraciones: Resultados: 
Preproceso: Condiciones de frontera, Gráficos, datos, 


Geometria, mallado monitoreos, condiciones convergencia, 
de resultados animaciones 


Ilustración 6. Resumen del proceso de CFD 
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METODOLOGIA 


El desarrollo del proyecto conlleva dos etapas: el disefio basado en principios 
hidráulicos y cálculos, así como la creación en un modelo 3D del disefio mediante un 
software de disefio asistido por computadora (CAD); y la simulación del modelo producido 
dentro de un ambiente de Fluido Dinâmica Computacional (CFD). 


Inicio 


Calcular la velocidad y 
diámetro específico para Proponer dimensiones del 
obtener un diámetro de impulsor y rendimientos 
disefio 


Obtener las 
condiciones de 
operación 


Calcular y disefiar Verificar Calcular desempefio 
la carcasa desempefio esperado del impulsor 


No 


Seleccionar material y 
configuración de 
simulación 


Realizar 
simulaciones 


Crear modelado 3D 
de la bomba 


Completar documentación Determinar Soluciones 
del diseio NSPHr adecuadas 


Ilustración 7. Diagrama de flujo de la metodologia propuesta. 


11 DISENO DE LA BOMBA 

Se considerará el disefio que pueda implementarse en un sistema de calentamiento 
de agua mediante energia solar para un centro deportivo. La bomba permitirá el flujo entre 
un árbol de calentadores solares que mantendrán caliente un tanque, del cual se podrá 
obtener el agua caliente para su uso. La mayor parte del agua caliente se utiliza en duchas, 
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cuyo consumo es de 20 L/min. Por lo tanto, un gasto pertinente sería de al menos 240 L/ 
min y una altura relativamente baja, pues solo se bombeará a través de los calentadores y 
otros componentes, se estima una altura de 20 m. 

Debido a los requerimientos, se puede acoplar la bomba a un motor eléctrico común, 
como el mostrado del catálogo de WEQ motores que permitirá una potencia nominal de 7.1 
kW a 2895 rom 


1.1 Cálculos iniciales para el disefio 


Lo primero que se debe especificar es la velocidad específica de la bomba, que 


utilizando una aproximación propuesta: 


nQ!2 


“ 53H3/4 e 


Os 
Posteriormente, se utiliza el diagrama de Cordier para obtener un diâmetro específico 
y, en consecuencia, obtener el diâmetro externo del impulsor. Quedando un diámetro inicial 
para el disefio de 14 cm, no obstante, se opta por un diâmetro más pequefio con la finalidad 
de disminuir el tamafio de la bomba para su aplicación. 
Marchegiani (2004) y Parra & Robles (2016) proponen medidas para el ángulo de 
salida de los álabes En la práctica, los ángulos más comunes para B, se encuentra entre 20 
y 35º. Ala entrada, B, tiene valores entre 50 a 15º. La velocidad meridiana en la salida se 
suele escoger por debajo de un 85 % de la velocidad meridiana en la entrada. En la tabla 1 
se especifican los cálculos iniciales de disefio. 


Parámetros de operación 


Eficiencias supuestas 


Caudal Q= 240 L/min 0.004 m?3/s — |Hidráulica Yn = 0.9 
Altura H = 20 m Volumétrica Mv = 0.96 
Potencia P 1100 W 825.2632 Mecánica Nm = 0.95 
Velocidad n 2895 rpm 48.25 rps Total 


D(gH)H* 
Ge E 


Tabla 1. Cálculos iniciales del disefio 


De acuerdo con la experiencia de Gómez y Zarate (2004), la transformación de la 
velocidad en altura de presión siempre está ligada a pérdidas de consideración. Por este 
motivo, con álabes curvados hacia delante es de esperar un rendimiento peor que con 
álabes curvados hacia atrás. 

A este punto es momento de seleccionar el número de alabes (Z), que de acuerdo 
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con la teoría de Pfleiderer se recomienda utilizar la expresión: 


14 


(% + 5) Bi + B> 
=D sen 


1.2 Obtención de las dimensiones del impulsor 


Para obtener las dimensiones, se utiliza un proceso de iteraciones donde se calcula 
la altura útil y el caudal producido a partir de una primera aproximación usando diâmetros, 
anchos y ángulos de los álabes a la entrada y salida, así como el número de álabes y su 
espesor basándose en los puntos previamente establecidos. 

Para determinar si un disefio de impulsor es adecuado, se deben calcular la altura 
útil y el caudal que generan, esto se logra determinando los triângulos de velocidades 
teóricos y haciendo los ajustes necesarios de eficiencia y otras correcciones. 

Primero se calculan las velocidades tangenciales con la ecuación (1) y, usando el 
supuesto de que es una bomba en condiciones óptimas a =90º se calculan las velocidades 
meridionales usando la ecuación (2), para posteriormente determinar el resto del triângulo 
con propiedades trigonométricas. Finalmente se utiliza la ecuación de Euler (4) para 
finalmente hacer las correcciones necesarias con una eficiencia hidráulica supuesta de 
n,=0.9 y utilizando el coeficiente de influencia de los álabes de la ecuación (9) y finalmente 
obtener una altura útil. De igual forma se calcula un caudal aproximado con la ecuación (2) 
y usando una eficiencia volumétrica n,=0.96 que es la más común para bombas radiales. 

Se recomienda el uso de una hoja de cálculo para permitir modificar los datos y 
conocer la forma en que cada uno impacta sobre el caudal y la altura. En la tabla 1 se 
muestran tres opciones de diseÃo que cumplen con la altura y el caudal deseado. 


1.3 Determinar la curva de los álabes 


Una forma sencilla de determinar la curvatura de los álabes es mediante un avance 
lineal de los mismos. Se segmenta en partes similares tanto el radio como el ángulo del 
álabe de la forma: 

t2 = 19 B> — B1 


Ars A =" 15 


n=n + Arj, B; = B; + AB) 16 


Utilizando estas ecuaciones se divide la curva en m partes enteras, y evaluando j 
con enteros entre O y m, se pueden obtener los valores correspondientes al ángulo B, que 
debe existir en el radio r, En la tabla se muestra el procedimiento con ayuda de una hoja 
de cálculo. 
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1.4 Disefio de la voluta 


Fernández Díez define la voluta de una bomba como un dispositivo que capta el 
líquido lanzado por el rodete impulsor y lo envía a la tubería de impulsión. Su cálculo y 
disefio se basan en el supuesto de que la componente circunferencial con de la velocidad 
en la câmara espiral, varía en proporción inversa al radio. 

Para determinar la forma de una voluta se considera que la velocidad varía 
inversamente al radio, se utiliza una espiral logarítmica para formar el contorno y se utiliza 
la sección transversal (sea rectangular o circular) para adecuar la manera en que la voluta 
recoge el fluido, haciéndolo pasar por una sección cada vez más amplia reduciendo su 
velocidad y, consecuentemente, aumentando la presión. 

Para esto, se define la circulación F como la velocidad a través del perímetro cerrado 
que abarca el rodete, que es una magnitud constante para una voluta dada en el régimen 
de trabajo. Por lo tanto, tenemos: 


CF =2nrc,=cte 17 


También se considera que el gasto q, en cualquier sección crece de manera 
proporcional al ángulo de inclinación O de la espiral que forma la voluta, desde la sección 
inicial que también es la sección final. 

0º 
Go = sen Qa 18 
360 


La teoría nos dice que, usando todas estas consideraciones, se puede obtener las 
dimensiones necesarias para el contorno de la voluta. Ya sea para una de forma rectangular 
con ancho constante b donde se requiere el radio del contorno; o una voluta con sección 


transversal circular con radio variable p, y con una distancia desde el eje del impulsor hasta 


el centro de la sección circular a. 


ERAS E pi 19 
360 210, (2) 


6º Qu, [280ar; 


= 20 
360 TJ 3607 


Pv = 
Variando el ángulo 6 desde O hasta 360 podemos obtener las dimensiones de la 
voluta, similar al modelo con el que se calculan las curvaturas de los álabes. Es importante 
notar que cuando 9=360º, r, toma el valor de R que es el radio máximo de la voluta, y en 
un perfil circular R=r,+2p, de igual forma con 0=360º, finalmente a=R-p, y se recomienda 
usar valores de r,. 
Fernández Díez (2012) y Khaing (2019) recomiendan que la relación - tome 
2 


valores entre 1.5 y 2 y la relación 5 entre 1.03 y 1.05. 
r2 
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21 SIMULACIÓN 


La simulación es la etapa final y más crítica del proceso, se considera el punto 
de inflexión de los resultados. A este punto debemos conocer suficiente sobre la teoría 
de funcionamiento de las bombas centrífugas para poder preparar la simulación con las 
condiciones óptimas, esto con la finalidad de facilitar la obtención de resultados, así como 
asegurar la calidad de los mismos 


2.1 Criterios para la simulación 


En orden para realizar una simulación adecuada es necesario establecer las 
condiciones de operación del sistema a simular, así como las condiciones de frontera. 

Es esencial conocer las propiedades del flujo, con una tubería de succión de alrededor 
los 42 mm de diámetro y con las propiedades del agua a 20ºC, se leen la viscosidad y 
densidad con 0.001 kg/m s y 1000 kg/m, respectivamente. Sabiendo que el volumen de 
trabajo es de 0.004 mº/s, podemos confirmar que el número de Reynolds establece un flujo 
turbulento. Como dato adicional para la simulación, la rugosidad común del hierro fundido 
es de 120 micras. 

En cuanto a condiciones de operación, lo primero que debe considerarse es que la 
bomba trabaja con un flujo volumétrica constante de 0.004 mº/s el cual se mantendrá a la 
entrada y salida en tuberías de succión y descarga. Esto se puede corroborar a partir de un 
equilibrio de masa; el flujo másico a la succión y a la descarga deben ser la misma, esto se 
transcribe a que el flujo métrico en respectivos puntos también es lo mismo puesto que la 
densidad del fluido de trabajo cambia extremadamente poco como para ser considerado. 

Las bombas trabajan con una sobrepresión a la descarga con la finalidad de evitar 
la cavitación, además de facilitar la entrada del fluido al impulsor donde se le imparte 
realmente la energía. Por lo tanto, se considera una sobrepresión de 5 metros de columna 
de agua, que se traduce como una presión absoluta de 150375 Pa. 
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n= 2895 rpm 303.1637 rad/s 


Propiedades del agua 
n= 0.001 kg/m s 


Q = 0004 m'3/s 
Vo = 2.887165 m/s 


Re = 121260.9 
Pom = 49050 Pa 5m.ca 
. Po = 150375 Pa 


Tabla 2. Resumen de los datos para simulación 


Finalmente se considera la operación del impulsor con una velocidad de rotación 
de 2895 rpm, lo cual se logra simulando una rotación dentro del espacio del impulsor 
y manteniendo fijo el impulsor. A su vez se debe excluir la voluta de esta rotación 
reconociéndolo como un estator, es decir una parte que no gira. 

A continuación, se resumen los puntos antes mencionados: 


1. Un volumen de entrada de 0.004 mº/s en la tubería de succión. 
2. Una presión de entrada de 150375 Pa en la tubería de succión. 
3. Una rotación del fluido en el impulsor de 2895 rpm. 

4. Paredes reales de estator en la voluta. 


5. Un caudal de salida en la tubería de succión de 0.004 m%/s 


2.2 Simulación en Solidworks 


El primer modelo sobre el que se realizará la simulación es en Solidworks, que 
es un software de disefio asistido por computadora comúnmente utilizados por alumnos 
de ingeniería mecánica pues ofrece un ambiente amigable y relativamente sencillo para 
resolver problemas de disefio, simulación y es a su vez compatible con otros sistemas. 

Para realizar la simulación en este programa, basta con activar el aditamento de 
Flow Simulation incluido y abrirlo dentro del archivo donde se quiere efectuar el estudio. 
La forma más sencilla de realizarlo es mediante un asistente que permite configurar 
rápidamente el proyecto. Este asistente comenzará por definir el nombre y detalles del 
proyecto, y procederá a solicitar la información necesaria, que incluye el tipo de simulación 
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(interna o externa, estacionario o transitorio, etc.), dominio del fluido, los materiales y 


sus propiedades. Después de Ilenar la información, se procederá con la introducción de 


condiciones de frontera. 


Wizard - Project Name 


= É Input Data 
E Computational Domain 
ME Component Control 
CB Fluid Subdomains 
Ef Boundary Conditions 
[as Fans 
MD Heat Sources 
“9 Porous Media 
SG Initial Conditions 
= Goals 
“QB Local Initial Meshes 
= Dg Results 
[HE Mesh 
2% Cut Plots 
< Surface Plots 
E Isosurfaces 
= Flow Trajectories 


Ilustración 8. Asistente para proyectos de Flow Simulation de SolidWorks 


Project 
Project name Projectft] 


Comments: 


Configuration to add the project 


Configuration Use Current 


Configuration neme: Bomba centrifuga 


Cancel 


Help 


Antes de continuar con el uso de la herramienta, es crucial el incluir otros elementos 


al diseio que facilitarán determinar las condiciones de operación y mejorarán los resultados. 


En este caso se crean dos tubos que simulan la entrada y la salida del fluido, esto para 


obtener espacios donde se puedan aplicar las condiciones de la entrada del fluido, así 


como monitorear los resultados a la salida. También se agregó un cuerpo que comprende el 


dominio que estará girando constantemente, un sólido que comprende el espacio existente 


entre las paredes del impulsor, con la finalidad de simular el movimiento del fluido en esa 


región. 
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a 


Ilustración 9. Cuerpo que simula la rotación (blanco) 


Finalizado el modelo, el asistente permitirá configurar las condiciones de frontera 
y de operación. Para este momento es importante contar con todas la información 
previamente establecida. Se comienza con el tipo de estudio que se quiere hacer, si a 
través del tiempo o estacionario, posteriormente se Ilenan las propiedades del fluido, así 
como verificar la geometria del sistema, creando cavidades donde habrá flujo y cerrando el 
dominio computacional. 

Oportunamente los sistemas de simulación cuentan con bibliotecas previamente 
cargadas con las propiedades de los materiales, así como una extensa documentación. 
En caso de SolidWorks esto facilita el inicio de la simulación. Sin embargo, existen campos 
que el usuario debe conocer para ser rellenados a consciencia, como lo son la rugosidad 
de las paredes. 

Una vez terminada la configuración se deben insertar las condiciones de frontera y 
operación. Las condiciones en la frontera son las que presenta el sistema por el universo 
que lo rodea, provienen externamente, en este caso la presión a la entrada. De operación 
son las que dictan como se comportaran el resto de las variables en la configuración, para 
lo estudiado se considera la velocidad de giro del impulsor y el flujo volumétrico del punto 
de operación que se quiere estudiar. 
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Además, existen otras condiciones que son importantes, como las paredes que 
existen y como estas se mueven. Por ejemplo, las paredes de la voluta son estáticas 
(estatores) y en contraste las paredes del impulsor giran. Estas características afectan los 
resultados al momento de considerar la viscosidad del fluido. 

En conjunto a las condiciones que se conocen de la simulación, también se deben 
establecer objetivos o zonas donde se medirán las características del fluido. Debido a que el 
caso de estudio es una bomba centrífuga, se monitorearán las presiones y las velocidades 
a la salida, esto permitirá hacer cálculos para el desempefio de la bomba. Teniendo de esta 
forma la presión de entrada y de salida, se puede obtener la altura útil producida por el 
modelo. Todas las configuraciones que se agreguen al proyecto aparecerán en el panel de 
navegación al costado izquierdo de la ventana de SolidWorks. 


£5; bombal 
E & Input Data 
ds | Computational Domain 
fed BB Fluid Subdomains 
- a & Rotating Regions 
8 Rotating Region 1 
x mi Boundary Conditions 
E a] Inlet Volume Flow 1 
did “Bm Outlet Volume Flow 1 
fee “Bm Environment Pressure 1 
| du “ Real Wall 1 
E a Goals 
ato A SG Average Static Pressure 1 
dera é BA SG Average Total Pressure 1 
fes o SG Average Dynamic Pressure 1 
ca EB SG Volume Flow Rate 1 
loga é SG Average Velocity 1 
na 8, SG Average Velocity (V) 1 


8 Mesh 


pia né Global Mesh 
E)-- :-Bg Results (1.fld, 5.000 s) 


Ilustración 10. Panel de navegación del proyecto 


Otra característica y principal para obtener los mejores resultados posibles es el 


mallado. Para los usuarios inexpertos, Solidworks ofrece realizar un mallado automático en 
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caso de no ser previamente definido. El mallado divide el dominio en sectores poligonales 
y cada sector se estudiará como un elemento. Al final el espacio se habrá dividido en un 
número finito de elementos a los cuales se les aplicarán las ecuaciones que determinan las 
propiedades del flujo en cada elemento y posteriormente las del siguiente elemento y así 
consecuentemente. 

Mediante un proceso iterativo, un algoritmo de soluciones aplica las condiciones 
previamente determinadas al mallado y cada uno de los elementos de este calculando 
paso por paso las características y reportando constantemente hacia los resultados que se 
desean monitorear. Este proceso se repite hasta que el algoritmo converge a un conjunto 
de resultados numéricos, el cual da la solución y permite obtener los resultados y la 
información que el usuario desee. 


Time = 4.000 s 


Ilustración 11. Mallado sencillo de SW 
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Outlet Volume Flow 
0.004 mº3is 


AngularVelocity | 
-303.163691 radis 


InletVolume Flow | 
0.004 m3is 


ço < 


Environment Pressure | 


150375 Pa 


— 


Ilustración 12. Diagrama isométrico de la simulación en SW. Incluye: 


X 


Volumen de entrada (Inlet volumen flow) 
Volumen de salida (Outlet volumen flow) 
Velocidad angular (Angular velocity) 
Pared real de estator (Real Wall Stator) 


Presión de entrada (Environment Pressure) 


2.3 Simulación en ANSYS Fluent 


La simulación mediante ANSYS Fluent es similar a Solidworks si las comparamos 
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desde un punto de vista meramente estructural. Sin embargo, existen diferentes 
metodologias y configuraciones más complejas que pueden implementarse dentro de la 
interfaz de Fluent. 

Lo primero que resalta es que Fluent maneja distintos programas para cada uno 
de los pasos principales para realizar una simulación, todo esto se puede acceder desde 
la plataforma principal, Workbench. Tiene un programa de dibujo 3D, otro para generar el 
mallado, el tercero es un asistente que permite configurar las condiciones iniciales, cuenta 
con una consola que soluciona el problema y finalmente una interfaz para presentar y 
exportar los resultados. El producto de un programa es exportado al siguiente; por ejemplo, 
la geometria se exporta para realizar el mallado sobre ella. 


w A 


E és, Fluid Flow (Fluent) 


2 |] Geometry 4 
3 É ven 4 
4 (8 setup A à 
5 Solution 4 
6 O results A 4 


7 tpa Parameters pa 
Fluid Flow (Fluent) 


Ilustración 13. Etapas de la simulación en Fluent. 


SpaceClaim es el programa de dibujo 3D integrado a Fluent. Permite crear sólidos, 
superfícies, cuerpos revolucionados, planos y demás geometria que posteriormente puede 
ser nombrada para utilizarse en la siguiente etapa. 

Meshing toma las geometrias y facilita la configuración de las condiciones 
importantes para generar el mallado, como lo son el número de celdas, la geometria de 
estas, el mallado dinámico, las zonas críticas, si se trata de un fluido o sólido y si ciertas 
geometrias deben excluirse del mallado. 
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0.000 0.050 0.100 (m) x 
E = | 


0.025 0.075 


Ilustración 14. Mallado del modelo 3 en Meshing 


Fluent es el componente principal, que cuenta con asistentes para determinar 
las condiciones de frontera, propiedades del fluido, método de solución, condiciones de 
solución, variables de control, las propiedades a monitorear; así como inicializar la solución. 
También, dentro de la misma plataforma se calculan las soluciones del modelo y las 
presenta activamente dentro de la interfaz hasta encontrar los resultados de acuerdo con la 
configuración propuesta por el usuario. 

CFD-Post extrae los resultados de Fluent para ser observados y presentados de 
distintas formas, ya sea como hoja de cálculo, gráficos, reportes. También permite exportar 
la información obtenida a otros formatos, así como generar capturas y cálculos posteriores 


a la simulación. 
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llustración 15. Condiciones para la simulación en Fluent: 
Azul: Entradas de presión manométrica de 49050 Pa y flujo de 0.004 m3/s 
Amarillo: Paredes del impulsor girando a 2895 rpm 
Blanco: Paredes estáticas 


Rojo: Salida de flujo y presión atmosférica 


202082 
ACADEMIC 
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ANÁLISIS DE RESULTADOS 
11 PRESENTACIÓN DE LOS RESULTADOS TEÓRICOS DE DISENO 


1.1 Resultados téoricos de las dimensiones del impulsor 


En la tabla 2 se pueden observar las iteraciones de disefo, englobando las 
dimensiones del impulsor, el triângulo de velocidades y el desempefio esperado de la 
bomba. En cuanto al radio de entrada del impulsor se tiene que r,= 0.021 m y el radio de 
salida r,= 0.070 m. En el ángulo de álabe para la entrada es de 8 = 20º y en la salida de 8,= 
10º, y por último se especifica el ancho del impulsor, en la entrada se tiene que b,= 0.020 
my en la salida b,= 0.018 m. Para los álabes en esta iteración se tiene que el numero Z= 3, 
con un espesor e= 0.005 m, con t = 0.65384 y t,= 0.7955. Se resume dentro de esta tabla, 
los resultados de los triângulos de velocidades en la entrada y salida. Para la velocidad 
absoluta de entrada se tiene que c,= 2.3172 m/s y para la salida de c,=17.6323 m/s. En la 
velocidad tangencial de entrada se tiene el valor de u = 6.3664 m/s y en la salida u,= 21.22 
m/s. Para los resultados de la velocidad relativa se tiene que en la entrada el valor de w = 
6.7750 m/s, y en la salida w,= 13.3442 m/s. Para la velocidad de la proyección tangencial 
relativa en la entras, se tiene que w, = 6.3664 m/s y en la salida w, = 3.6006 m/s. En la 
proyección tangencial en la velocidad absoluta de la entrada, c, = 0.0 m/s y en la salida de 
Coy= 17.62 m/s. En cuanto a la velocidad meridional en la entrada se tiene que c, = 2.3172 
m/s y en la salida de c, = 0.6349 m/s. Para los triângulos de velocidades se tiene que los 
ângulos son los siguientes: 8,= 20º, a,= 90º, y,= 70º, B,= 10º, a=2.06º y y,= 167.94º. 

El desempefio esperado para la bomba (ver tabla 2), se resume a los siguientes 
puntos: altura de Euler (H )= 38.11825 m, coeficiente de Píleiderer (p= 0.950343), altura 
corregida (H )= 22.84555 m, caudal final (Q)= 0.004165 mº/s, altura final (H)= 20.561 m y 
potencia final (P)= 921.1103 W. 

Para el primer disefio, se utilizaron las dimensiones del impulsor como las 
establecidas en la primera interación de disefio. Para el segundo diseio se efectuaron 
cambios en el radio de entrada (rl= 25 mm) y salida (r,= 70 mm), en el ângulo de entrada 
(B,= 25º) y de salida (8,= 15º), así como en el ancho de entrada (b = 12 mm) y de salida 
(b,= 10 mm), se efectuaron también cambios en el número de álabes, incrementando 2 
más, teniendo un total de Z=5 (ver tabla 2). Para el tercer disefio se efectuaron cambios en 
el radio de entrada (r1= 22 mm) y salida (r,= 70 mm), en el ángulo de entrada (8 = 20º) y 
de salida (8,= 15º), así como en el ancho de entrada (b = 20 mm) y de salida (b,= 15 mm), 
se efectuaron también cambios en el número de álabes, incrementando 1 más, teniendo un 


total de Z= 4, respecto a la primera iteración (ver tabla 2). 
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Primera iteración de disefio 


Dimensiones del impulsor 
Dimension | mtrada | sanida | Atabes 
Radio E 
Ángulo álabe 0.005 m 
Ancho . = ! 0.65384 Tz = 0.7955 
Velocidad Entrada Ángulos 
Absoluta =— 2.3172 m/s = 17.6323 m/5 20º 
Tangencial = 6.3664 m/s 21.22 m/s 
Relativa 6.7750 m/s 13.3442 m/s 
Proyección tang| Wiu = 6.3664 m/s = 3.6006 m/s 
Proyección tang| Ciu = 0.00 m/s == 17.62 m/s 
Meridional 2.3172 m/s 0.6349 m/s 
Desempefio esperado de la bomba 

Altura de euler H,= 3811825m Caudal Final =" 0.004165 m/3/s 

Coef. Pfleiderer PV = o0.950343 Altura Final 20.561 m 

Altura corregida Hz = 22.84555m Potencia final 921.1103 W 


Primer disefio 
Dimensiones del impulsor 


Dimensión | Entrada | salida | Álabes 


70 mm |Número Z= 
10 * Espesor = 


18 mm 


Segundo disefio 
Dimensiones del impulsor 


Dimensión | entrada | *  salida | Álabes 


70 mm |Número Z= 5 
15 * Espesor e = Sm 
10 mm 


Tercer disefio 


Dimensiones del impulsor 

Dimensión | Entrada | **Salida | Álabes 
i 70 mm |Número 7Z= 
15 * Espesor e = 

15 mm 


Tabla 3. Proceso de obtención de dimensiones y dimensiones finales 


1.2 Resultados teóricos de las curvas de álabes 


Para determinar las curvas de álabes se tomaron como referencia las configuraciones 
de los 3 disefios propuestos en tabla 2. En la tabla 3, se resumen los resultados obtenidos 
para los 3 disefios. En el diseÃo 1 se tiene que ri= 15 mm, r2= 100 mm, Ar= 17 mm, B,= 
25º, B,= 10º y AB = -3. Para los álabes del disefio 2 se tiene que que ri= 15 mm, r2= 95 
mm, Ar= 16 mm, B,= 35º, B,= 15º y AB= -4. Para los álabes del disefo 3 se tiene que rt= 
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20 mm, r2= 100 mm, Ar= 16 mm, B,= 30º, B= 15º y AB=-3. 


Curvas de álabes 


Tabla 4. Tabla curva de álabes 


1.3 Resultados teóricos del disefio de la voluta 


Los resultados de la voluta se muestran en la tabla 4. En la tabla se resumen las 
propiedades de la voluta, que incluyen la proyección tangencialc, = 17.39 m/s, la circulación 
[= 7.65 m?/s y el caudal impulsor Q = 0.00384 mº/s. En cuanto a las dimensiones de la 
voluta se tiene que el radio impulsor r2= 70 mm, el radio máximo R= 85.09 mm, el ancho de 
la voluta b= 25.00 mm y el radio inicial r3= 75 mm. 


DISENO DE LA VOLUTA 
EEE 
Do jo [xo 


Dimensiones de la voluta | 90 [003155) 7740 
70 mm 135 


85.09 mm 


75 mm 


ss 
360 


Tabla 5. Disefio de voluta 
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21 PRESENTACIÓN DE LOS RESULTADOS DE DISENO 


Los disefios presentados a continuación, se realizaron en el software Solidworkse 
a partir de los parámetros teóricos calculados en el apartado anterior. La ilustración 16, 
muestra el modelo tridimensional en isométrico de la voluta, donde se tomaron como 
referencia las dimensiones de la voluta y la curva de la voluta, especificados en la tabla 
4. La ilustración 17, muestra los disefios de los tres impulsores con base a los resultados 
teóricos de la tabla 3 y 4. En la ilustración 18 se muestra el plano de la voluta con las 
dimensiones especificas del disefio. Las ilustraciones 19, 20 y 21, corresponden a los planos 
de los 3 impulsores disefiados. Se anexan los planos especificados anteriormente, como 
parte de posibles réplicas de los componentes y su ensamblaje, para fines de simulación y 
obtención de otros parámetros no contemplados en este trabajo de investigación, o en su 
defecto mejoras en el disefo. Los planos de talles se elaboraron a través de Solidworksê 
Drawing Tools. En la ilustración 22 se muestran los ensambles del impulsor y la voluta de 
los disefos 1,2 y 3. 


A 


Ilustración 16. Modelo tridimensional de la voluta 
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A. 


ad Es 


Ilustración 17. Modelos tridimensionales del impulsor. De izquierda a derecha; disefio 1,2 y 3 
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SECCIÓN F-F 


Titulo 
VOLUTA 
Fecha: Dibuló: 
| riam camoménico, camosécror 
e 


Ilustración 18. Plano de la voluta 
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SECCIÓN D-D 


SECCIÓN E-E 


MODELO IMPULSOR 1 


Fecha: Dibujó: 
FEB 2021 CASTRO HUÉRIGO, CARLOS HÉCTOR 


Ilustración 19. Plano del modelo del impulsor 1 
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Nie 


. 
est 


Y 

77) 8 
ZA 
4) 


CNN 


SECCIÓN D-D 


SECCIÓN E-E 


MODELO IMPULSOR 2 


Escala: Fecna: Jó: 
12 | FEB2021 CASTRO HUÉRIGO, CARLOS HÉCTOR 


Ilustración 20. Plano del modelo del impulsor 2 
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SECCIÓN D-D 


SECCIÓN E-E 


MODELO IMPULSOR 3 


Escda: Fecho: Dibujá: 
1:2  FEB2021] CASTRO HUÉRIGO, CARLOS HÉCTOR 


llustración 21. Plano del modelo del impulsor 3 
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Ilustración 22. Ensamblaje de la voluta y los impulsores; disefio 1,2 y 3 


31 PRESENTACIÓN DE LOS RESULTADOS DE SIMULACIÓN 


Se midieron las propiedades de salida a la descarga, la presión total, dinâmica, 
estática, la velocidad vertical y la velocidad absoluta. Los resultados se contrastan a 
continuación: 

Se hace énfasis en la presión estática a la salida, pues esta puede ser considerada 
como la presión que un manómetro mediría a la descarga, de esta forma podemos obtener 
la altura producida por la bomba como la diferencia de estas presiones en metros de 
columna de fluido. Para la simulación en Solidworks, la presión medida es absoluta, es 
decir, considera la presión atmosférica; en Fluent, es manométrica y no considera la presión 


de atmósfera. Por este motivo se tiene que realizar una corrección a la presión estática y 
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total. Dicha información se contrasta con los resultados esperados de los cálculos teóricos 
y se presentan en gráficos comparativos. 

En conjunción, se analizarán gráficas de distribución de las presiones y velocidades 
dentro de cada uno de los modelos. Mediante un análisis detenido del comportamiento del 
fluido en un gráfico del alzado de la bomba, se pueden detectar anomalías, puntos críticos, 
oportunidades de mejora del disefio. En la tabla 6 se observan los resultados de presión 
estática, presión total, presión dinámica, flujo volumétrico, velocidad absoluta y velocidad 
en Y, obtenidos de SolidworksO Flowsimulationy ANSYSO Fluent, para los disefios 1,2 y 8. 
En la ilustración 23 se muestran los gráficos de presión absoluta obtenidos de Solidworks 
Flowsimulationy ANS YSO Fluent, como parte de la comparativa entre ambos softwares. En 
la ilustración 24 se observa el gráfico de contraste de alturas producidas por la bomba en los 
disefios 1,2 y 3, en donde se diferencian los resultados obtenidos por ambos simuladores. 
Las ilustraciones de 25 a 30, muestran los resultados de simulación de velocidad y presión 
para los modelos 1, 2 y 3, analizados con Solidworksê Flowsimulation. Las ilustraciones de 
31 a 36, muestran los resultados de simulación de velocidad y presión para los modelos 1, 
2y83, analizados con ANSYSO Fluent. 


RESULTADOS DE SIMULACIÓN 


V | 
H Propiedad Unidades 


EE E E 
Presióntotal [| 4277775/ moszas|pa) | 
restóninânico | sen nel, 


Fujovolumétdico | o volumétrico | -0004 (o) | -0004 004 Ima3/s] | s 
Velocidad absoluta | 2.187239 2.1813| [m/s] 
Velocidad en Y 2.056307 1.4277/[m/s] 


3912179] 295230/[pa) | 


Disefio 1 


Presióntotal | 3038243] 297751 B)jpa) | 


Presión dinámica 2606.20: 2521.8/[Pa] 


Flujo volumétrico 


Velocidad absoluta | 2.28334 
Velocidad en Y 2.05630 


Presión total 378638.2| 314123.8 O] 
| 


Tabla 6. Resultados de Simulación 


Disefio 2 


Disefio 3 
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CONTRASTE DE ALTURAS 


| Diseio [ 1 | 2 | 3 [unidades 
[20.56 | 21.73 | 22.58 | 


20.56 | 2173 | 22.58 
sw| 28.06 | 24.58 | 23.02 
Fluent| 26.36 25.12 26.81 m 


Tabla 7. Contraste de resultados entre modelos de simulación y teoria 


GRÁFICA 
| É RE 
Estática | 2 [39121785 396555 | [Pal | 


absoluta [3 [37598073 413085 | [Pal | 
ontota | o2rrrrs asas | (Pal | 
Presondel [2 [s9382406 3950768 | (Ps) | 
| 3 [37382 aisamss | (pa) | 


Tabla 8. Presiones absolutas de simulación 


Presión Estática Absoluta 
430000 


420000 


410000 


400000 
390000 
380000 
370000 
360000 
350000 
1 2 3 


m=SW =Fluent 
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Presión Total Absoluta 


440000 
430000 
420000 
410000 
400000 
390000 
380000 
370000 
360000 
350000 
1 2 3 
m=SW s=Fluent 
Ilustración 23.Gráficos de contraste Presiones Absolutas 
Altura Total [m] 
30.00 


25.00 

20.00 

15.00 

10.00 

5.00 

0.00 
1 2 3 


mTeoria mSW =Fluent 


llustración 24. Gráfico de contraste altura total producida por la bomba 


Observando cada uno de los disefios se puede dictaminar una cercanía entre la 
simulación que proporciona Solidworks y la de Fluent. Según ilustración 23, el disefio 3 
parece ser el de mayor diferencia entre los valores de presión, realmente esto se traduce 
en aproximadamente 3 m de altura final según gráfico de ilustración 24. Esta pérdida de 
presión y ganancia de altura se documentó también en Gómez & Zarate (2004). se puede 
observar que en los modelos 1 y 2, las diferencias de presión son menores (ver ilustración 
283), resultando una discrepancia de 2 m de altura final de dichos modelos (ver ilustración 24). 
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Es importante resaltar que los resultados teóricos, fueron los que arrojaron los 
resultados más bajos en cuanto la ganancia de presión por unidad de masa. Este fenómeno 
ocurre a menudo a nivel teórico, ya que cuando el modelo es más conservador, ofrece una 
forma más sencilla y rápida de realizar las estimaciones, por lo que es de esperarse que 
los resultados no sean tan cercanos a la realidad, por lo tanto, bajo esta metodologia, los 
sistemas y máquinas son disefiados para sobrepasar las características deseadas. 

En cuanto a los resultados en conjunto de cada uno de los disefios se mantuvieron 
alrededor de los 24 m, sobrepasando los 20 m que se esperaban al realizar las 
aproximaciones iniciales, validando de esta manera los 3 disefios analizados. 


344609.61 
320554.28 
296498.95 
272443.62 
248388.30 
224332.97 
200277.64 
176222.31 
152166.98 
128111.65 
104056.32 
80000.99 
55945.66 
31890,34 
7835.01 


Relative Pressure [Pa] 


Cut Plot 1: contours 


L. 


Ilustración 25. Simulación SW. Modelo 1, gráfico de presiones 
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21.219 
19.704 
18.188 
16.672 
15157 
13.641 
12.125 
10.610 
9.094 
7.578 
6.063 
4.547 
3.031 
1.516 
0 


Velocity [mis] 


Cut Plot 1: contours 
Cut Plot 2: contours 


La 


Ilustración 26. Simulación SW. Modelo 1. Gráfico de velocidad 
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469299.44 
431380.50 
393461.55 
355542.61 
317623.67 
279704.72 
241785.78 
203866.84 
165947.90 
128028.95 
90110.01 

52191,07 

1427212 

-23646.82 
-61565.76 


Relative Pressure [Pa] 


Cut Plot 1: contours 


L.. 


Ilustración 27. Simulación SW. Modelo 2. Gráfico de presiones 
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21.782 
20.227 
18.671 
17415 
15.559 
14.003 
12.447 
10.891 
9.335 
7779 
6.224 
4.668 
3112 
1.556 
Li) 


Velocity [m/s] 


Cut Plot 1: contours 
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Ilustración 28. Simulación SW. Modelo 2. Gráfico de velocidad 
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300984.34 
279361.49 
257738.64 
236115.79 
214492.94 
192870.09 
171247.23 
149624.38 
128001.53 
106378.68 
84755.83 
63132.98 
41510.13 
19887.28 
-1735.57 


Relative Pressure [Pa] 


Cut Plot 1: contours 


La 


Ilustración 29. Simulación SW. Modelo 3. Gráfico de presiones 
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velocity [mis] 


Cut Plot 1: contours 
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Ilustración 30. Simulación SW. Modelo 3. Gráfico de velocidad 
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Pressure 
presiones 


3.053e+05 
2.875e+05 
2.698e+05 
2.521e+05 
2.343e+05 
2.166e+05 
1.988e+05 
1.811e+05 
1.634e+05 
1.456e+05 
| 1.279e+05 
+ 1.102e+05 
| 9.242e+04 
» 7.468e+04 
' 5.695e+04 
É 3.921e+04 
* 2.147e+04 
[Pa] 


Ilustración 30. Simulación Fluent. Modelo 1. Gráfico de presiones 
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Ilustración 31. Simulación Fluent. Modelo 1. Gráfico de velocidad 
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Pressure 
presiones 


2.914e+05 
2.738e+05 
2.561e+05 
2.384e+05 
2.207e+05 
2.031e+05 
1.854e+05 
1.677e+05 
1.501e+05 
1.324e+05 


[ 1.147e+05 
+ 9.704e+04 


7.937e+04 
6.170e+04 


| 4.403e+04 
| 2.635e+04 
- 8.683e+03 


[Pal 


Ilustración 32. Simulación Fluent. Modelo 2. Gráfico de presiones 
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Velocity 

velocidades Figure 2 
2113e+01 
1.981e+01 


1.849e+01 
1.717e+01 


“ 


Ilustración 33. Simulación Fluent. Modelo 2. Gráfico de velocidad 
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Pressure 
presiones Figure 2 
3.080e+05 
2.900e+05 
2.720e+05 
2.541e+05 
2.361e+05 
2.181e+05 
2.001e+05 
1.821e+05 
1.641e+05 
1.461e+05 
+ 1.281e+05 
+ 1.101e+05 
| 9.211e+04 
1 7412e+04 
: 5.612e+04 
É 3.813e+04 
FÉ 2014e+04 


[Pa] 


Ilustración 34. Simulación Fluent. Modelo 3. Gráfico de presiones 
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Velocity 

velocidades e dr 
2.040e+01 
1.912e+01 
1.785e+01 
1.657e+01 
1.530e+01 
1.402e+01 
1.275e+01 
1.147e+01 
1.020e+01 


im s*-1) | 


Ilustración 35. Simulación Fluent. Modelo 3. Gráfico de velocidad 


Observando las gráficas presentadas se pueden notar que los fenómenos resaltados 
en el marco teórico se presentan también en las simulaciones. Como se comentó, los 
gráficos permiten, de forma llamativa y amigable, discernir lo que está sucediendo dentro 
del modelo simulado. 

Si seguimos el flujo del líquido en los resultados de simulación, lo primero que 
notaremos es que las presiones más bajas se encuentran a la entrada del impulsor. Esto es 
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un fenómeno esperado pues es en esta área donde el fluido tiene altas velocidades y por lo 
tanto sufre de altas pérdidas por fricción y conversión a energia cinética derivado del efecto 
Bernoulli, documentado también por Oliveira (2002), Evans (2012), Lal & Deshmurkh (2018), 
Gúner & Ózbayer (2020) y Altyntsev et al. (2021) . En estas zonas se debe monitorear los 
niveles de presión para asegurar que la presión no baje a más del punto de vaporización y 
de esta forma se presente la cavitación. 

También podemos analizar como la velocidad absoluta dentro del impulsor aumenta 
respecto al radio desde donde se mide. De esta forma corroboramos el funcionamiento del 
impulsor, ya que este debe transmitirle cada vez más energía al fluido conforme pasa más 
tiempo siendo impulsado por este (Chapallaz et al., 1992). 

Otro aspecto importante que se puede notar en la cuchilla presente entre el cono 
difusor y el punto donde comienzo la voluta. Este espacio puede causar turbulencias y ser 
fuente de pérdidas, siempre ha sido un punto importante de diseho (Duarte, 2011). 

La voluta es una sección que no se considera tan relevante en comparación al 
impulsor, sin embargo, es un componente de gran impacto en la eficiencia de la bomba 
como sistema. En los gráficos se puede ver cómo, a pesar de que la voluta cada vez recoge 
más material alrededor de la trayectoria que la espiral define la carcasa, la velocidad no 
tiene variaciones importantes, lo cual contribuye a menor cantidad de pérdidas y por lo 
tanto mayor eficiencia. 

Siguiendo la misma idea de la voluta, el cono difusor convertirá la velocidad recogida 
por la carcasa en presión. En la tubería de descarga es el punto donde la presión debe ser 
mayor, ya que, para ese punto, la energia ya fue suministrada por el impulsor y también fue 
convertida a presión por la carcasa. Dicha zona es un buen indicador de la eficiencia de la 
bomba, pues a ese punto se pueden medir las características del fluido. Efectivamente, y 
como se dijo anteriormente, midiendo la entrada y la salida se puede obtener una altura útil. 
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CONCLUSIONES 


El disefio de una máquina es un problema que puede ser abordado desde distintos 
puntos de vista, como puede ser uno meramente práctico y empírico, o como un ejercicio 
de cálculos complejos que se aproximan a la realidad. En este proyecto se demostró 
la viabilidad de una metodologia que toma las fortalezas de dos métodos de disefio 
fundamentales; la teoría y la simulación. 

La facilidad de aplicar la bibliografia respecto a disefio de bombas centrífugas, así 
como su expresión rápida de los resultados mediante hojas de cálculo, son puntos fuertes 
que tienen ventaja sobre las consideraciones teóricas; no obstante, estos procedimientos 
llevan consigo muchas suposiciones, valores aproximados y tienden a desviarse de la 
realidad. La Dinámica de Fluidos Computacional (CFD) ofrece un entorno más cercano a la 
realidad, con cálculos más complejos pero especializados dentro de modelos complicados 
y de muchas variables. Su desventaja principal es que requieren más recursos tanto 
computacionales como de tiempo de cálculo. 

Al realizar una primera aproximación mediante cálculos teóricos con la bibliografia 
y suposiciones adecuadas, podemos reducir la cantidad de modelos simulados necesarios 
para obtener un disefio óptimo para las condiciones de operación requeridas. 

Se demostró que los resultados de ambos métodos son cercanos entre sí, esto 
fortalece la idea de una metodologia unificada. Finalmente, los gráficos, animaciones y 
demás recursos que ofrece la CFD, son eficaces para detectar problemas en los disefios, 
y así analizar problemas que puedan surgir y corregirse antes de elaborarse un prototipo 
físico. 

En cuanto al disefio y simulación de los modelos podemos resumir lo siguiente: 

* | Dentro del disefio 1 se presentan pérdidas de presión por la cuchilla de la voluta. 


* | Enlatubería de descarga del modelo 2, la presión generada no fue tan grande, 
en comparación a las pérdidas por la cuchilla de la voluta. 


* | Enel modelo 3 la voluta convierte gradualmente la velocidad en presión de for- 
ma correcta, al igual que el cono difusor. La excepción es el espacio ubicado en 
la transición voluta-cono, donde parece haber pérdidas de presión. 


* | Los modelos 1 y 2 las presiones bajan significativamente a la entrada, por lo 
que se recomendaria aumentar la altura disponible. 


* | Eldisefo 3 parece tener el mejor desempefio de los que se presentaron 


*- Fluent permitió dar resultados más exactos y con mayor precisión que Soli- 
dWorks. 


* La descarga en los tres modelos genera turbulencias que disminuyen el rendi- 
miento de la bomba. 
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RECOMENDACIONES 


El proyecto realizado permitió dar un acercamiento más integral de los diversos 
modelos utilizados dentro del disefio de una máquina compleja, en este caso de estudio 
se propuso una bomba centrífuga. Comenzando por un primer acercamiento con cálculos 
y bibliografia que se recomienda ampliamente en un modelo metodológico, iterativo y con 
variables de disefio a elegir para encontrar una inicialización cercana a la realidad del 
disefio a producir. Posteriormente se canalizó la información obtenida a un entorno virtual 
de simulación de dinámica de fluido computacional que permite acercar aún más a la 
realidad las características deseadas del modelo. 

La metodologia, las ecuaciones, la bibliografia, así como las hojas de cálculo y 
ambos softwares de simulación y disefio, se recomiendan ampliamente como herramientas 
fundamentales para el proceso de generar máquinas (o componentes de estas) capaces de 
trabajar a las condiciones de operación esperadas. 

Todo lo anterior mencionado puede corroborarse de acuerdo con la forma en 
que los modelos lograron acercarse entre sí con los resultados obtenidos. Si bien cada 
entorno tiene sus ventajas, por ejemplo, los cálculos teóricos son accesibles, no requieren 
de procesamiento de cómputo ni software y hardware costoso, sin embargo, los errores 
pueden ser más recurrentes en este tipo de técnicas. Por otra parte, tenemos métodos (con 
base a las condiciones teóricas) más precisos y complejos como ANSYS Fluent, sin olvidar 
un balance más amigable pero efectivo que representa SolidWorks, con la desventaja de 
que tenemos que tener acceso al software y hardware. De ahí la importancia de conocer 
ambos, con la intención de siempre tomar la mejor opción para resolver los problemas que, 
como ingenieros, afrontamos dentro de un ambiente de desarrollo tecnológico y práctico, 
pudiendo incluir soluciones combinadas de validación. 
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ANEXOS 


11 SW MODELO 1 


General Info 

Run at: LAPTOP-RL53F2F0 
Number of cores: 8 
Iterations: 4073 

Physical time: 5.000 s 

CPU time: 2036 s 

Log 


Mesh generation started O 


16:53:02, Nov 01 


Mesh generation normally finished O 


16:53:05, Nov 01 


Preparing data for calculation O 


16:53:05, Nov 01 


Calculation started O 


16:53:06, Nov 01 


Calculation has converged since the following 
criteria are satisfied: 4073 


17:27:08, Nov 01 


Max. phys. time is reached 4073 


Calculation finished 4073 


17:27:08, Nov 01 


Calculation Mesh 
Basic Mesh Dimensions 


Number of cells in X 8 
Number of cells in Y 16 
Number of cells in Z 10 
Number Of Cells 

Cells 7192 
Fluid cells 7192 
Irregular cells 0 
Trimmed cells 0 


Maximum refinement level: 2 
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Goals 


Name Unit Value Progress | Criteria Delta an 
SG Average 
iso Pressure Pa 425389.19 | 100 557981.69 18.6610352 On 
a i 427777.50 557981.598 17.413624 On 
SG Average 
Dynamic 2388.30 4.39000118 1.27658957 On 
Pressure 1 
E ad om -0.0040 1.30034807e-05 | 2.20817375e-12 | On 
o m/s 2.187 100 0.00211933035 | 0.000583320545 | On 
a m/s 2.056 100 3.65630235e-05 | 1.13517151e-09 | On 
Min/Max Table 

Name Minimum Maximum 

Density (Fluid) [kg/m/3] 997.55 997.56 

Pressure [Pal] 109160.01 445934.61 

Temperature [K] 293.20 293.25 

Temperature (Fluid) [K] 293.20 293.25 

Velocity [m/s] 0 21.219 

Velocity (X) [m/s] -21.054 20.977 

Velocity (Y) [m/s] -20.675 21.049 

Velocity (Z) [m/s] -4.184 2.179 

Velocity RRF [m/s] 0 17.757 

Velocity RRF (X) [m/s] -15.752 16.398 

Velocity RRF (Y) [m/s] -17.749 15.579 

Velocity RRF (Z) [m/s] -4,184 2.179 

Vorticity [1/s] 3.88e-07 2061.79 

Relative Pressure [Pa] 7835.01 344609.61 

Shear Stress [Pa] 0 171.38 

Botileneck Number [] 7.2393817e-10 1.0000000 

Heat Transfer Coefficient [W/m/2/K] 0 0 
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ShortCut Number [] 7.2183213e-12 1.0000000 
Surface Heat Flux [W/m/2] 0 0 

Surface Heat Flux (Convective) [W/m/2] -2.759e+09 3.609e+09 
Acoustic Power [W/m3] 0 2.933e-09 
Acoustic Power Level [dB] 0 34.67 


21 SW MODELO 2 


General Info 

Run at: LAPTOP-RL53F2F0 
Number of cores: 8 
Iterations: 4097 

Physical time: 5.000 s 

CPU time: 2299 s 

Log 


Mesh generation started O 


13:45:34, Nov 01 


Mesh generation normally finished O 


13:45:37, Nov 01 


Preparing data for calculation O 


13:45:37, Nov 01 


Calculation started O 


13:45:39, Nov 01 


Calculation has converged since the following 
criteria are satisfied: 4097 


14:23:57 , Nov 01 


Max. phys. time is reached 4097 


Calculation finished 4097 


14:23:58, Nov 01 


Calculation Mesh 
Basic Mesh Dimensions 


Number of cells in X 8 
Number of cells in Y 16 
Number of cells in Z 10 
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Number Of Cells 


Cells 6781 
Fluid cells 6781 
Irregular cells 0 
Trimmed cells 0 
Maximum refinement level: 2 
Goals 
; se Use in 
Name Unit Value Progress | Criteria Delta convergence 
SG Average 
Static Pressure | Pa 391217.85 | 100 2343084.87 49.6322255 On 
1 
SG Average 
Total Prosguród Pa 393824.06 | 100 2343085.03 51.5094807 On 
SG Average 
Dynamic Pa 2606.20 35 2.6556772 7.52576844 On 
Pressure 1 
ENA m3/s |-0.0040 | 100 1.26708927e-05 | 1.65423231e-12 | On 
Flow Rate 1 : i ' 
SG Average m/s 2.283 38 0.00116888345 | 0.00309669952 | On 
Velocity 1 
SG Average ; . 
Velocity (Y) 1 m/s 2.056 100 8.49353775e-05 | 8.50402859e-10 | On 
Min/Max Table 
Name Minimum Maximum 
Density (Fluid) [kg/m"3] 997.54 997.57 
Pressure [Pal] 39759.24 570624.44 
Temperature [K] 293.18 293.30 
Temperature (Fluid) [K] 293.18 293.30 
Velocity [m/s] 0 21.782 
Velocity (X) [m/s] -21.053 21.137 
Velocity (Y) [m/s] -20.792 21.049 
Velocity (Z) [m/s] -3.115 2.708 
Velocity RRF [m/s] 0 17.285 
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Velocity RRF (X) [m/s] -15.584 15.233 
Velocity RRF (Y) [m/s] -17.282 14.805 
Velocity RRF (Z) [m/s] -2.115 2.708 
Vorticity [1/s] 9.52e-05 20159.76 
Relative Pressure [Pa] -61565.76 469299.44 
Shear Stress [Pa] 0 1036.36 
Bottleneck Number [] 4.4787863e-11 1.0000000 
Heat Transfer Coefficient [W/m/2/K] 0 0 

ShortCut Number [] 7.5164547e-10 1.0000000 
Surface Heat Flux [W/m/2] 0 0 

Surface Heat Flux (Convective) [W/m/2] -2.938e+09 2.537e+09 
Acoustic Power [W/m3] 0 4.748e-09 
Acoustic Power Level [dB] 0 36.77 


31 SW MODELO 3 


General Info 

Run at: LAPTOP-RL53F2F0 
Number of cores: 8 
Iterations: 3728 

Physical time: 5.000 s 

CPU time: 1731 s 

Log 


Mesh generation started O 


10:42:37, Oct 30 


Mesh generation normally finished O 


10:42:41, Oct 30 


Preparing data for calculation O 


10:42:42, Oct 30 


Calculation started O 


10:42:43, Oct 30 


Calculation has converged since the following 
criteria are satisfied: 3728 


11:11:33, Oct 30 


Max. phys. time is reached 3728 


Calculation finished 3728 


11:11:34, Oct 30 


Warnings: A vortex crosses the pressure opening Boundary Condition : Environment 


Pressure 1; Inlet flow/outlet flow=24.6203 


Calculation Mesh 
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Basic Mesh Dimensions 


Number of cells in X 8 
Number of cells in Y 16 
Number of cells in Z 10 
Number Of Cells 
Cells 7083 
Fluid cells 7083 
Irregular cells 0 
Trimmed cells 0 
Maximum refinement level: 2 
Goals 
Ê FR Use in 
Name Unit Value Progress | Criteria Delta convergence 
SG Average 
Static Pressure Pa 375980.73 | 100 1327557.6 12.8979741 On 
1 
SG Average 
Total Pressure 1 Pa 378638.20 | 100 1327557.69 12.4916794 On 
SG Average 
Dynamic Pa 2657.47 100 10.4459428 1.90252327 On 
Pressure 1 
e Flow | mag/s | -0.0040 100 1.69727575e-05 | 1.75257638e-12 | On 
Ri PN GRADO m/s |2.305 100 0.00496553513 | 0.000812937647 | On 
Velocity 1 , : , 
SG Average m/s |2.056 100 8.18451026e-05 | 9.00958863e-10 | On 
Velocity (Y) 1 E ; , 
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Min/Max Table 


Name Minimum Maximum 
Density (Fluid) [kg/m"3] 997.55 997.56 
Pressure [Pa] 99589.43 402309.34 
Temperature [K] 293.20 293.24 
Temperature (Fluid) [K] 293.20 293.24 
Velocity [m/s] 0 21.219 
Velocity (X) [m/s] -21.053 20.977 
Velocity (Y) [m/s] -20.917 21.049 
Velocity (Z) [m/s] -5.193 1.941 
Velocity RRF [m/s] 0 15.799 
Velocity RRF (X) [m/s] -14.220 14.157 
Velocity RRF (Y) [m/s] -15.794 14.363 
Velocity RRF (Z) [m/s] -5.193 1.941 
Vorticity [1/s] 1.03e-05 1757.75 
Relative Pressure [Pa] -1735.57 300984.34 
Shear Stress [Pa] 0 909.82 
Bottleneck Number [] 1.8093594e-09 1.0000000 
Heat Transfer Coefficient [W/m/2/K] 0 0 
ShortCut Number [] 3.6067590e-11 1.0000000 
Surface Heat Flux [W/m/2] 0 0 

Surface Heat Flux (Convective) [W/m/2] -2.538e+09 6.150e+09 
Acoustic Power [W/m3] 0 5.392e-09 
Acoustic Power Level [dB] 0 37.32 
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41 ANSYS MODELO 1 


1. File Report 
Case FFF 
File Path | | CAUsersicarloiDocumentsYTHiResidenciasNTrabajoNFluenASimulacion fisimulaciont . 


filesidpONFFAFluenhFFF-1 


-00700.dat.gz 


File Date | 06 December 2020 
File Time | 01:47:32 PM 
File Type | FLUENT 

File Version | 20.2.0 


Table 1. File Information for FFF 


2. Mesh Report 


Domain Nodes Elements 
impulsor 11247 54389 
piezaí 7478 36799 
pieza2 1043 3640 
pieza 3 ensamblaje1 2192 1755 
All Domains 21960 96583 


Table 2. Mesh Information for FFF 


3. Physics Report 


Domain - impulsor 


Type 


solid 


Domain - piezaí 


Type 


cell 


Domain - pieza2 


Type 


cell 


Domain - pieza 3 ensamblaje1 


Type 


cell 


Table 3. Domain Physics for FFF 
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Domain Boundaries 
Boundary - contact region 2 contact region 3 contact region 4 contact . 
region 7 contact region 4 contact region 8 src impulsor 
Type INTERFACE 
impulsor 
Boundary - wall impulsor 
Type WALL 
Boundary - contact region 2 contact region 3 contact region 4 contact . 
region 7 contact region 4 contact region 8 src piezaí 
Type INTERFACE 
pizal Boundary - contact region 2 contact region 3 contact region 4 contact . 
region 7 contact region 4 contact region 8 trg piezaí 
Type INTERFACE 
Boundary - contact region 2 contact region 3 contact region 4 contact . 
region 7 contact region 4 contact region 8 src pieza2 
Type INTERFACE 
Boundary - contact region 2 contact region 3 contact region 4 contact . 
region 7 contact region 4 contact region 8 trg pieza2 
) Type INTERFACE 
pieza2 
Boundary - outlet flow 
Type MASS-FLOW-INLET 
Boundary - wall pieza2 
Type WALL 
Boundary - contact region 2 contact region 3 contact region 4 contact . 
region 7 contact region 4 contact region 8 trg pieza 3 ensamblajei 
Type INTERFACE 
pieza 3 Boundary - inlet flow 
ensaios Type MASS-FLOW-INLET 
Boundary - inlet p 
Type PRESSURE-INLET 


iter continuity x-velocity y-velocity z-velocity k omega velocidad- velocidad- presion-ab presion-di presion-es time/iter 
694 2.4827e-03 2.0460e-04 2.7525e-04 1.6044e-04 5.9202e-04 3.8372e-04 2.1899e+00 1.4277e+00 4.0887e+05 3.1711e+03 3.0754+05 0:00:05 7 


Reversed flow on 518 faces (79.4% area) of pressure-inler 14, 


695 2.6387e-03 2.1126e-04 2.8028e-04 1.6247e-04 5.8699e-04 3.9042e-04 2.1884e+00 1.4277e+00 4,0884405 3.1680e+03 3.07528+05 0:00:05 5 
Reversed flow on 527 faces (82.0% area) of pressure-inlet 14. 

696 2.7055e-03 2.0995e-04 2.7918e-04 1.6948e-04 5.9890e-04 3,9370e-04 2.1869e+00 1.4277e+00 4.0885e+05 3.1650+03 3.07538+05 0:00:04 4 
Reversed flow on 516 faces (75.2% area) of pressure-inler 14. 

697 3.2105e-03 2.2886e-04 3.0212e-04 1.6830e-04 6.0075e-04 3.8970e-04 2.1855e+00 1.4277e+00 4.0879e+05 3.1621e+03 3.07468+05 0:00:03 3 


Reversed flow on 515 faces (80.0% area) of pressure-inlet 1 
€S8 2.57808-03 2.0467e-04 2.6984e-04 1.71308-04 5 


uy 


e-04 3.8959e-04 2.1841e+00 1.42778+00 4.08782+05 3.15928+03 3.07458+05 0:00:02 


Reversed flow on 516 faces (79.7% area) of pressure-inlet 14. 


699 2.3937e-03 2.0147e-04 2.6220e-04 1.6366e-04 5.99]5e-04 3.8842e-04 2.18272+00 1.4277e+00 4.08728+05 3.1563e+03 3.0739e+05 0:00:01 1 
Reversed flow on 521 faces (81.1% area) of pressure-inlet 14. 
700 2.3303e-03 1.9431]e-04 2.6750e-04 1.6033e-04 5.9331e-04 3.8720e-04 2.18138+00 1.42778+00 4.08708+05 3.1535e+03 3.07378+05 0:00:00 O 


Registering ReportDefFiles, ("C:lUsersicarlolDocumentsNITHlResidenciasTrabajolFluentNSimulacion 1llsimulacionl filesldp0NFFFFluent,.Ypresion-din-rfile.out 
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5|I ANSYS MODELO 2 


1. File Report 
Case FFF 
File Path CAUsersicarloDocumentsNTHiResidenciasNTrabajoNFluenASimulacion ASimulacion2 | 
filesidpONFFAFluentFFF-2-00500.dat.gz 
File Date | 06 December 2020 
File Time | 05:43:23 PM 
File Type | FLUENT 
File Version | 20.2.0 


Table 1. File Information for FFF 


2. Mesh Report 


Domain Nodes | Elements | Tetrahedra | Wedges | Pyramids | Hexahedra | Polyhedra 
impulsor 6856 32669 32669 0 0 0 0 
pieza1 4596 21499 21499 0 0 0 0 
pieza2 1046 3666 3666 0 0 0 0 
REI aq | 2192 1755 0 60 0 1695 0 
All Domains | 14690 59589 57834 60 0 1695 0 
Table 2. Mesh Information for FFF 
3. Physics Report 
Domain - impulsor 
Type solid 
Domain - piezaí 
Type cell 
Domain - pieza2 
Type cell 
Domain - pieza 3 ensamblaje1 
Type cell 
Table 3. Domain Physics for FFF 
Anexos 
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Domain Boundaries 


Boundary - contact region 2 contact region 3 contact region 6 src 


impulsor 
Type | INTERFACE 
impulsor 
Boundary - wall impulsor 
Type | WALL 
Boundary - contact region 2 contact region 3 contact region 6 trg piezaí 
Type | INTERFACE 
piezal 
Boundary - contact region 7 src 
Type | INTERFACE 
Boundary - contact region 2 contact region 3 contact region 6 src pieza2 
Type INTERFACE 
Boundary - contact region 2 contact region 3 contact region 6 trg pieza2 
Type | INTERFACE 
pieza2 
Boundary - outlet flow 
Type | MASS-FLOW-INLET 
Boundary - wall pieza2 
Type | WALL 
Boundary - contact region 7 trg 
Type | INTERFACE 
pieza 3 Boundary - inlet flow 
ensamblaje! Type | MASS-FLOW-INLET 
Boundary - inlet p 
Type PRESSURE-INLET 
Reversed flow on 1 face (0.2% area) of pressure-inlet 14. 
Reversed flow on 1 face (0.2% area) of pressure-inlet 14. 
Reversed flow on 1 face (0.2% area) of pressure-inlet 14. 


500 2.1737e-03 6.5216e-05 1.3235e-04 6.2459e-05 3.0045e-04 1.6741e-04 1.4277e+00 


1.7990e+00 2.52188+03 2.95238+05 3.96552+05 0:00:00 [o 
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61 ANSYS MODELO 3 
1. File Report 


Case 


FFF 


CAUsersicarloDocumentsNTHiResidenciasNTrabajoNFluenhSimulacion 3isimulacion 


File Path | 5 filesidpONFFAFluenhFFF.1-5-02500.dat.gz 
File Date | 02 December 2020 

File Time | 08:55:00 PM 

File Type | FLUENT 


File Version | 20.2.0 


Table 1. File Information for FFF 


2. Mesh Report 


Domain Nodes Elements 
impulsor 12513 62320 
piezal 10012 48614 
pieza2 1041 3626 
pieza 3 ensamblaje1 2192 1755 
All Domains 25758 116315 


Table 2. Mesh Information for FFF 


3. Physics Report 


Domain - impulsor 


Type 


solid 


Domain - piezaí 


Type 


cell 


Domain - pieza2 


Type 


cell 


Domain - pieza 3 ensamblaje1 


Type 


cell 


Table 3. Domain Physics for FFF 
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iter continuity x-velocity 
2495 1.9143e-03 6.10278-05 
2496 1.8 

2497 
2458 
2499 


presion-ab ve velocidad- presion-di presion-es presion-to 
4.1309e+05 00 1.7803e+00 2.36326+03 3,1176e+05 3.1413+05 
2.3634e+03 


635e+03 
6378+03 
638e+0 


Note: zone-surface: cannot create surface from sliding interface zone. 


Table 4. Boundary Physics for FFF 


Domain Boundaries 
Boundary - contact region 2 contact region 3 contact region 8 src 
impulsor 
Type INTERFACE 
; Boundary - contact region 4 contact region 9 src impulsor 
impulsor 
Type INTERFACE 
Boundary - wall impulsor 
Type WALL 
Boundary - contact region 2 contact region 3 contact region 8 trg piezaí 
Type INTERFACE 
Boundary - contact region 4 contact region 9 src piezaí 
piezal 
Type INTERFACE 
Boundary - wall impulsor contact region 2 contact region 3 contact . 
region 8 trg 
Type INTERFACE 
Boundary - contact region 2 contact region 3 contact region 8 src pieza?2 
Type INTERFACE 
Boundary - contact region 2 contact region 3 contact region 8 trg pieza2 
Type INTERFACE 
Boundary - inlet ps 
pieza2 
Type WALL 
Boundary - outlet flow 
Type MASS-FLOW-INLET 
Boundary - wall pieza2 
Type WALL 
Boundary - contact region 4 contact region 9trg 
Type INTERFACE 
pieza 3 Boundary - inlet flow 
ensamblaje 1 Type MASS-FLOW-INLET 
Boundary - inlet pe 
Type PRESSURE-INLET 


time/iter 
0:00:13 


5 
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